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DIE DYNAMISCHE THEORIE DER GEZEITEN AUF 
EINEM MACLAURIN’SCHEN ELLIPSOID 

VON 

PKOF. D”- KARL HILI.LBRAND 


VORGELEGT IN DER SITZUNG AM 19. DEZEMBER 1912 


Über die dynamische Theorie der ozeanischen Gezeiten sind in den letzten Dezennien sehr ein¬ 
gehende fundamentale Arbeiten erschienen. Es sind dies insbesondere die umfangreichen Arbeiten 
G. H, Darwin’s, die in den »Proceedings of the Royal Society« und in der »Encyclopaedia Britannica« ver¬ 
öffentlicht wurden, ferner S. S. Hough’s: »On the Application of Harmonie Analysis to the Dynamical 
Theorie of the Tides«, I & 11, Lond. Phil. Trans. A. vol. 189 und 191, und die einschlägigen Kapitel in 
Poincare’s »Le 9 ons de mecanique celeste«, t. IIl (Theorie des Marees). 

Allen diesen Untersuchungen liegt die Annahme einer sphärischen Gleichgewichtsfigur zugrunde, 
eine Annahme, die streng genommen im Widerspruch steht mit der das Phänomen wesentlich mitbestim¬ 
menden Rotationsgeschwindigkeit. Durch die Berücksichtigung der entsprechenden strengen Gleich- 
gewichtsfigur wird allerdings in den grofBen Zügen der theoretischen Darstellung wenig geändert und sie 
erscheint vielleicht umso weniger von wesentlicher Bedeutung, als ja doch der tatsächliche Verlauf der 
Gezeiten in überwiegender Weise durch die spezielle Konfiguration der Meeresbecken mitbestimmt wird, so 
zwar, daß er mit dem aus gewissen einfachen Voraussetzungen gewonnenen theoretischen Verlauf nur in 
den allgemeinsten Umrissen übereinstimmt. 

Es muß aber doch die — auch in dieser Hinsicht — konsequente 'Einbeziehung der Rotations¬ 
geschwindigkeit als eine immerhin einmal durchzuführende Ergänzung dieser Theorie um so mehr ange¬ 
sehen werden, als gewisse Komponenten der Gezeitenbewegung sich dabei etwas anders darstellen, 
insbesonder das eine Laplace’sche Theorem eine Modifikation erfährt. Außerdem ist die Annahme einer 
der Rotation entsprechenden Gleichgewichtsfigur, also im aktuellen Falle eines Maclaurin’schen Ellipsoides 
durchaus kein Umstand, der die harmonische Analysis in irgendeiner Weise schwieriger gestaltet, so daß 
eigentlich kein Anlaß vorliegt, von einer mechanisch völlig korrekten Annahme abzugehen. Es würde 
sogar bei einer noch allgemeineren Annahme — der einer ellipsoidischen Gleichgewichtsfigur überhaupt 
— der gleiche analytische Gang mit Zuhilfenahme der Lame’schen Funktionen möglich sein. 

Die folgende Untersuchung soll sich auf den einzig aktuellen Fall eines Rotationsellipsoides 
beschränken. 

Zugleich mit dieser Verallgemeinerung soll eine möglichst übersichtliche Darstellung der dynami¬ 
schen Theorie der Gezeiten überhaupt gegeben werden und aus diesem Grunde auf numerische Details, 
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soweit sie nicht die erwähnte Ivorrektur betreffen, verzichtet und bezüglich eingehenderer quantitativer 
Angaben auf die oben zitierten Fundamcntalwerke verwiesen Averden. 

Noch soll Erwähnung finden, daß bei der Behandlung der beiden Hauptgieichungen wesentlich der 
von S. S. Hough eingeschlagene Weg verfolgt wurde und die Auflösungsmethode eine Verallgemeinerung 
der von ihm gegebenen ist. 


I. Die Differentialgleichungen der Gezeitenbewegung. 

Die rechtwinkligen Koordinaten eines Massenelementes der Flüssigkeitsschicht seien 5, C, bezogen 
auf ein System, das unveränderlich mit dem festen Teil des Himmelskörpers verbunden ist. Dieser rotiere 
mit einer Winkelgeschwindigkeit oj, deren Komponenten nach den rechtwinkligen Axen p. r sein mögen. 
Die Geschwindigkeitskomponenten des Flüssigkeitsteilchens sind dann 


Vx — — 

Vy = C 

= tj+pr^~qi 

wo S', r/, die Geschwindigkeiten relativ zum starren Körper sind. 

Ebenso hat man für die Komponenten der Beschleunigung 

Ux = Vx-^qVx — r Vy 

iiy — Vy+r Vx—p Vz 
= vi-^p Vy-~qVx. 

Man kann für die folgenden Untersuchungen annehmen, daß die Rotationskomponenten konstante 
Größen sind, das heißt also v'on Rotationsstörungen und etwaiger Abweichung der instantanen Drehungs¬ 
achse von der kleinsten Hauptträgheitsachse absehen. 

Dann erhält man, wenn für den Augenblick 

^ .2 + ^2 ^ 

gesetzt wird 

11X = {q C^~r r/) -f- co p p cos (to, p) —of- i 

Uy — 2 {r^—p'^) -h CO q [j cos (co, p)~ co- */] 

1 /.^ = C^'-h 2 (p r/ ^q^') -h CO r p cos (co, p) — co- C- 

usw. sind die Beschleunigungen der relativen Bewegung, 2 (q 'Cj — rr{) usw. sind die Komponenten 
der Coriolis’schen Kraft und ( 0 /; p cos (co, p)—oj-$ usw. die negativen Komponenten der Zentrifugalkraft. 
Das Potential der letzteren sei mit Q bezeichnet. 

Auf das Flüssigkeitsteilchen wirkt die Attraktion der Gesamtmasse des Himmelskörpers: des festen 
Teiles vermehrt um die der flüssigen Masse in der momentanen Konfiguration, das Potential derselben sei 
F; außerdem sei noch eine außerhalb gelegene Masse die als ausdehnungslos angenommen Avird, 
vorhanden, deren Attraktion als störende Kraft mit dem Potential 5 dazutritt. Das Gesamtpotential ist dann 

F= ß-DF+5 

so daß Q-hF=f/ das Potential der Schwere im betrachteten Punkt und 

v= s 

ist. 


i 
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Setzt man die Flüssigkeit als inkompressibel und reibungslos voraus und ihre Dichte gleich der 
Einheit, so sind die hydrodynamischen Gleichungen: 

r' +2(<?C'-rr,0= -^{V-p) 

ox 

und analoge für und p ist der hydrostatische Druck im betrachteten Punkt. 

Im Problem der Gezeiten gibt man diesen Gleichungen eine etwas andere Form, die dadurch bedingt 
ist, daß die tatsächlichen Lagen der Punkte sich nur wenig von den Gleichgewichtslagen unterscheiden 
werden, das heißt jenen Lagen, die sie nur unter dem Einfluß von U und ohne Eigenschwingung annehmen 
würden. 

Bezeichnet man die Koordinaten der Gleichgewichtslage mit z, die Abweichungen infolge der 
Gezeitenbewegung mit fv, so daß 


^ = X -h n, ri = y -y V, C = - + iv, 

so lauten nun die Gleichungen 


d- ii ^ ( d w 
+ 2 - 




dt 


d- V , / d n 
- \r 2 [r - 


df^ 

d'w 


+ 2 \p 

df^ \ dt 


dt 

dv 


r . 

dt j 

d IV 

p - 

dt 

d n 
dt 


= (r^ lA 

öx 


— iv^p) 

dy 


= (V P)- 


Dazu kommen noch die Kontinuitätsgleichung und die Grenzbedingungen. 

Bei der weiteren Behandlung dieser Gleichungen kann man sich sehr weitgehende Vereinfachungen 
erlauben, die beim Problem der Ebbe und Flut einer bedeckenden Flüssigkeitsschichte durchaus gerecht¬ 
fertigt sind. 

Vor allem kann man die Verschiebungskomponenten ?/, v, iv, um deren Ermittlung es sich zunächst 
handelt, als kleine Größen betrachten, von denen nur erste Potenzen berücksichtigt werden sollen. 

Daraus folgt in ei'ster Linie, daß die Euler’sche und Lagrange’sche Betrachtungsweise identisch 
werden, weil sich beide um Größen zweiter Ordnung unterscheiden, wie aus 

d n' 0 0 /F , 0 , 0 /F , 

_ . — -}_- ^ - y ^ - ly 

dt 0/ 0.r dy dz 


unmittelbar hervorgeht. 

Ferner können in der Störungsfunktion, die ja von derselben Ordnung wie die Deviationen ist, 
unmittelbar die Koordinaten der Gleichgewichtslagen eingeführt werden. 

. Was das Potential der Schwere U für ein Massenelement der deformierten Flüssigkeitsschichte 
anbelangt, so kann dasselbe mit der gleichen Annäherung durch folgende Betrachtung gefunden werden. 
Es sei Uq sein Wert auf der ungestörten Oberfläche, deren Gleichung demnach 


= Konst. 

ist. Die vermöge der Gezeiten gestörte Oberfläche wird von dieser nur durch Abstände getrennt sein, die 
erster Ordnung sind. Denkt man sich von den Randpunkten eines Plächenelementes der Gleichgewichts¬ 
oberfläche Normalen bis zum Durchschnitt mit der wirklichen Oberfläche gezogen, so kann also dieses 
kleine Prisma nach der obigen Annahme wie ein Volumelement behandelt und die ganze Deformation 
kann durch das Hinzufügen dieser kleinen Massen gedacht werden, die positiv oder negativ zu nehmen 
sind, je nachdem die Normalen nach außen oder innen gerichtet sind. Ist C die Höhe dieses kleinen 
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Prismas, so kann — als Ort seines Schwerpunktes und damit auch als Ort dieses Deformationselementes 


angesehen werden. Das Potential der Schwerkraft in einem solchen unterscheidet sich nun in zweierlei 
Weise von 17^: erstens geht es, wenn man nur die Verschiebung C dieses einen Elementes in Betracht 
zieht, offenbar in 


Uq + 



oC 


über, wo die Ableitung die nach der Richtung der Normalen von bedeutet; dann tritt aber noch in dem 
Teil P von U die durch die Deformation verursachte Änderung der Attraktion, das ist also die Gesamt¬ 
attraktion der Deformationselemente — unter Berücksichtigung ihres Vorzeichens — dazu. Dieser Bestand¬ 
teil von P ist demnach gegeben durch 


rfJ. 


WO äo ein Obeiflächenelement, A seine Entfernung vom betrachteten Flüssigkeitsteilchen und das Integral 
über die ganze GleichgewichtsoberOäche zu nehmen ist. 


Bemerkt man weiter, daß ~ nichts anderes als die Acceleration g der Schwere an der 

l3«/o 

betreffenden Stelle der Fläche ist, aber eine Konstante, daher für die Potentialfunktion bedeutungslos 
ist, so hat man 






C ist die Strecke der Normalen der Gleichgewichtsoberfläche zwischen dieser und der freien Ober¬ 
fläche, genügt daher der sich auf die gesamte Deformation beziehenden Kontinuitätsgleichung 


jC^/o:=0 


das Integral über die ganze Oberfläche genommen. 

In den Bewegungsgleichungen hat man also zu setzen 

F— + It- + S. 

Die Unbekannten u, v, iv sind Funktionen der Zeit und der Koordinaten x,y, z. Den ersteren Zusam¬ 
menhang werden die ßewegungsgleichungen ergeben, den letzteren die Kontinuitätsgleichung und die 
Grenzbedingungen. Die Art des funktionellen Zusammenhanges mit der Zeit ist aber sofort angebbar. 

Die linken Seiten der Gleichungen sind lineare Funktionen der Ableitungen nach der Zeit mit 
konstanten, nur von x,y,z abhängigen Koeffizienten, die rechtsstehende Funktion U~p-\-S ist im ersten 
Teil U~p eine lineare Funktion von n, v, iv gemäß der hier festgesetzten Annäherung, mit ebensolchen 
Koeffizienten, S ist eine bekannte Funktion der x,y, z und der Zeit, welch letztere Abhängigkeit man als 
eine Exponentialreihe voraussetzen kann. Derartige Gleichungssysteme werden aber bekanntlich durch 
ebensolche Exponentialreihen integriert, so zwar, daß man ein partikuläres Lösungssystem in der Form 
voraussetzen kann: 

11 = V = IV = - 

Die Uq, Vq, Wq sind Funktionen der Gleichgewichtslage, also von x,y, z; X ist entweder identisch mit 
dem Faktor t in einem Glied der Reihe S (erzwungene Schwingung) oder ergibt sich in ganz bestimmter 
Weise aus dem von S unabhängigen Teil des Integrals (freie Schwingung). Der letztere enthält die will¬ 
kürlichen Konstanten. 
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Für ein derartiges partikuläres Integral hat man also 

du . 

dt de- 

Setzt man F—/> = X“ tL, so müssen die gesuchten Größen als Funktion des Ortes sich durch die 
Ableitungen dieser Funktion ^ in der Weise bestimmen, daß 

Ihi + 2\iqw-rv) — X^ 

dx 

0 

X-v + 2X (ru — piv) — — 

‘dy 

)ew -y 2X {pv —qn) — X- . 

02: 

Es ist also das Problem auf die Bestimmung einer einzigen Funktion ^ reduziert, die aus anderen 
Bedingungsgleichungen gewonnen werden muß. Es sind dies die Kontinuitätsgleichung 

0?^ 0r' "dw 

— H-!-=0 

0.V 0^' 02: 

und gewisse Grenzbedingungen: daß für die freie Oberfläche 

p — Konst. 

und für die festen Teile der Begrenzung die zu ihr normale Verschiebungskomponente Null ist. 


2. Die Gleichungen für Flüssigkeitsschichten geringer Tiefe. 

Dieser Fall, der einzig und allein bei den ozeanischen Gezeiten auf unserer Erde in Frage kommt, 
läßt eine gewisse Vereinfachung der hier auftretenden Relationen zu. Es soll angenommen werden, daß 
die im allgemeinen veränderliche Tiefe k der Flüssigkeitsschicht von der Größenordnung der Ver¬ 
schiebungen ist. 

Das an sich ganz willkürliche Koordinatens 3 ^stem werde nun so gewählt, daß die ;rjV-Ebene parallel 
zur Tangentialebene an dem betrachteten Punkt der Gleichgewichtsoberfläche ist. Für die als eben anzu¬ 
nehmende Umgebung des Punktes auf der Fläche wird die kleine Größe h eine gewisse Funktion von .r 
undjk sein. (Bei der in Wirklichkeit ganz willkürlichen Veränderlichkeit von h kann es sich natürlich nur 
um den Verlauf gewisser Mittelwerte handeln.) Wenn nun h überall eine kleine Größe ist, so gilt dasselbe 

auch für das Fortschreiten dieses mittleren Wertes nach irgendeiner Richtung, das heißt für und 

0.r 0^^ 

Berücksichtigt man die Grenzbedingung für den Meeresgrund, vermöge welcher die normale VerschiebungS' 
komponente Null sein muß, also 

07f 0// 

«0 --1- -b = 0, 

dx 0j' 


so folgt daraus, daß die Komponente für die untere Begrenzung eine Größe zweiter Ordnung ist. Da 
die Änderung von w bis zur freien Oberfläche 


der Voraussetzung gemäß eine Größe derselben Ordnung ist, so kann die Komponente w überhaupt ver¬ 
nachlässigt werden, die dritte Bewegungsgleichung entfällt und die beiden anderen enthalten nur mehr die 
Rotationskomponente 


V — (0 cos 6, 
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wo 6 der Winkel der Normalen mit der Rotationsachse ist. Die Bewegungsgleichungen lauten dann: 

cr^ n 
dr^ 


dv 0 

2 oj cos 6 — = 

dt 0 .r 


-h 2 CO cos Ö- —-X- 

dr^ dt 


0 'p 
0y 


und man hat, dem obigen Resultat entsprechend^ 

)s2 


n =: 


+ 4 co‘‘^ cos 6 


0']> 2o3 0^M 

—— + cos 6 —- 

0.V X 0jr/ 


X2 

V = - 

X'*^ + 4 cos 6 



2 to - 0 fL 

-cos Ö —- 

X 0.V 


In dieser Lösungsform ist also bereits eine der Grenzbedingungen enthalten. 

Zur weiteren Behandlung der Aufgabe sollen für die Gleichgewichtsoberfläche krummlinige 
Gauß’sche Koordinaten |x und v eingeführt werden, so daß die Gleichung dieser Fläche in der Form 
gegeben ist 

^ =/l (|J-, v) y=f.i (p., v) S =/s ([>., v), 

und zwar sollen die Größen zwei Scharen orthogonaler Trajektorien vorstellen, wonach ein Linienelement 
ds der Oberfläche ausgedrückt ist durch 

ds = s/ Ed^^ + Gdv"^^ 

wo E und G die bekannten Gauß’schen Größen sind. 

Alle in den Gleichungen vorkommenden abhängigen Veränderlichen können nun als P'unktion dieser 
Koordinaten gedacht werden. 

Die Komponenten und u seien nun auch nach den durch |x und v angegebenen Richtungen orien¬ 
tiert, dann ist 


X-.i + 4 (ü-^ cob^ Ö 

X2 + 4ö,2cos2 6 \\/G0v 


+ — COS 6 — 

X G 0 V y 

2(0 0 6 \ 

— — cosö— 7 ^=!— 


Die Kontinuitätsbedingung kann auf Grund der Kleinheit der Tiefe in folgender Weise formuliert 
werden: Es sei h die Höhe der ungestörten Flüssigkeitsschicht an einem bestimmten Punkt der Ober¬ 
fläche. Betrachtet man den Raum eines unendlich dünnen rechtwinkeligen Prismas, dessen Grundlinien 
v/ifjjx, \/~Gdv und die entsprechenden gegenüberliegenden Linienelemente sind und dessen Höhe h ist, 
so ist durch die Seitenfläche \/Ed\Gi in einer bestimmten Zeit die Flüssigkeitsmenge uh\/'Ed]^ durch¬ 
gegangen und durch die gegenüberliegende 

0 

V h </E J [Ji -h (vh\/ Ed (x) d v, 

^ 0 V 

es resultiert also der Zuwachs 

—- {vJi\/E) r/[xr/v 

0 V 


und ebenso von dem darauf normalen Seitenflächenpaar 

— —^ {u hs/G) dv d |x. 
0 [X 
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Diese V^ermehrung des P'lüssigkeitsquantums drückt sich aber in der Größe 'Q, der Erhebung der 
freien Oberfläche über die ungestörte, aus und ist gegeben durch \/ÄGC^(xr/v; es ist daher 

— inlt\/G)+ — (vlis/E) = v/£GC 

8 [X 8 V 

und das ist die Kontinuitätsgleichung für den Fall kleiner Tiefen. Führt man die Darstellung der Defor¬ 
mationen n und V durch die Funktion ^ ein und setzt 




-4- 4 to^ cos^ 9 

so lautet die Kontinuitätsbedingung 


7 fi 7 

: cos 6 — h^ 

X 


\/EGi: = 


8(x 


K 


G ^ 
E 8 (X 


+ 


8 
8 V 


'V G 


8 

T7 


oder 


Dazu kommt noch die Bedingung der freien Oberfläche 

p ~ Konst. 

Xdo 




^^2 > 


+ 


8/7. 


+ S —Konst. 


8r{; 
8 V 


dk^ ^ 8 ^ 
8v 8p. 


oder, da für ein Potential additive Konstante irrelevant sind, 







Für S ist der der Größe X entsprechende Term Ce^^ zu substituieren für die erzwung^enen Schwin¬ 
gungen — C ist eine Funktion von [x und v und den Bewegungselementen des störenden Körpers — oder 
es ist S =: 0 und X bezieht sich auf die freien Schwingungen. Im letzteren Falle ergibt sich X aus einer 
bestimmten Bedingungsgleichung. 

Aus der Kontinuitätsbedingung und der Bedingung der freien Oberfläche sind die einzelnen partiku¬ 
lären Integrale für ^ und C abzuleitcn. Diese Größen sind — abgesehen von dem von der Zeit abhängigen 
Faktor—Funktionen des Ortes auf der Gleichgewichtsoberfläche, also gewisse Funktionen von (x und v, und 
die Lösung wird sich durch den Umstand besonders übersichtlich gestalten, daß sich derartige Funktionen 
— wenigstens für die hier in Frage kommenden Gleichgewichtsfiguren — in Reihen gewisser Elementar¬ 
funktionen entwickeln lassen, deren Integraleigenschaften eine sehr einfache Bestimmung der Koeffizienten 
zur Folge haben. Für die Kugel sind dies die Kugelfunktionen und für Ellipsoide die Lame’schen 
Funktionen. 

Da für das vorliegende Problem nur Maclaurin’sche Ellipsoide in Betracht kommen können, so 
sollen die weiteren Schritte zur Lösung desselben unter Zugrundelegung dieser speziellen Gleichgewichts- 
figur durchgeführt werden. 


3. Das Maclaurin’sche Ellipsoid. 

Die hier vorausgesetzte Gleichgewichtsform sei also ein Rotationsellipsoid. Die Lage eines Punktes 
auf der Oberfläche sei durch die drei elliptischen Koordinaten 

P> T 

gegeben, die im Falle eines Umdrehungskörpers auftreten, so daß der Punkt durch den Schnitt dreier sich 
orthogonal schneidender Flächen gegeben ist: des Rotationsellipsoides ^ 
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wo c die lineare Exzentrizität der Meridianellipse ist, des einschaligen Rotationshyperboloides, dessen 
Meridianschnitte mit dieser Ellipse konfokale Hyperbeln sind: 


-h y- 






= 1 


und der durch den Winkel cp bestimmten Meridianebene. Führt man statt die elliptische Kugelkoordinate ^ 
ein durch 

— “ sin xr, 
c 

wo also ö’ der halbe Öffnungswdnkel des zum H^^perboloid gehörigen Asymptotenkegels ist, so sind die 
rechtwinkeligen Koordinaten dieses Punktes 

p . • ft • 

X — - Sin 'S = p sin ^ sin 's 


y — —^ cos = p sin ü- cos cp 


1 


2: =1 - \yp- — C-\/c-- lX- — \/p' — C^ cos = p cos S cos \)‘, 

c 


wenn — = sin s, das heißt, s der Exzentrizitätswinkel ist 
P 

Daraus ergibt sich für das Linienelement ds der Ausdruck 

ds- — —^— (1—sin^ £ sin^^) dp^ + p^ (1—sin^ s sin^ ö-) 4- p- sin^ ^Jcp-, 

cos'^ £ 

also zerlegt nach den Komponenten in der Normalenrichtung und den beiden nach dem Meridian orien¬ 
tierten in der Tangentialebene. 

Ist ferner 6 der Winkel der Normalen mit der :3c-Achse, der hier mit dem Komplement der Polhöhe 
identifiziert werden kann, und r die geozentrische Distanz, so ist 


tg 6 cos £ tg ^ und r — psy 1 — sin- £ cos‘^ 0-. 

Soll das Ellipsoid bei der Winkelgeschwindigkeit w eine Gleichgewichtsfigur sein, so muß 

3 k-M £ (3 —2 sin-£) —3 sin £ cos £ 

10^ — —-. -- 


2 


sin^ £ 


sein, wo k die Konstante der Gravitation und M die Masse des Himmelskörpers bedeutet. 

Aus der diesem Werte entsprechenden Zentrifugalkraft und den bekannten Formeln für die Altraktions- 
komponenten der Rotationsellipsoide ergibt sich als Beschleunigung der Schwere auf der Oberfläche in 
einer Breite, die der Koordinate d' entspricht, 




3 M 


■ = —-^ (tg S-- s) -^ \/1 — sin- £ sin^ il 

p- sin^£ 


oder, wenn man die Beschleunigung am Pol mit öO bezeichnet, 

g = go\/l—sin^ s sin^ 9-. 

W'as nun die oben erwähnten zur Lösung heranzuziehenden Entwicklungsfunktionen anbelangt, so 
sind diese für elliptische Koordinaten die Lame’schen Funktionen, die aber für Rotationsellipsoide in die 
bekannten Kugelfunktionen übergehen. 
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Die Lame’schen Funktionen genügen der Gleichung 

if) = 0, 

wo den Laplace’schen Operator 

02 02 02 

— - + - Hh - 

d 8j|/^ 8 

bedeutet, transformiert für das betreffende Koordinatensystem: p, p, v. Ist dieses orthogonal, so zwar, daß 
das Linienelement ds gegeben ist durch 


ds"^ — a? d p d]y^ dv"^ 


so ist der Operator bekanntlich 


^2 _ 


aßY 


ßY 3 


8 / 7 a 8 


8p\a 8py 8(jL\ß 8[x/ 8v 




aß 8 


Im Falle der allgemeinen elliptischen Koordinaten ist ein Punkt gegeben durch drei Flächen zweiten 
Grades 




2 




p2 — ^2 


— 




= 1 




_yö _y-s 

WO p jJi> Z? > V > 0 ist. Die rechtwinkeligen Koordinaten sind dann 


— ^ y — v/p" - n/p-" - ^ s/c^-v 


, y _ 

hc b\^d^ — h‘^ 

Setzt man für den Moment 




= (p2_||,2) ([,2-v2) ^p2-vä) 

— (p2-c2) 

52 = (p,2_&2) (c2-[l,2) 

O — (&2_v2) (C2_v2) 


SO ist 


a 


ß = 


Q 


^\/p."-v2 

Q 

5\/p2-v2 


,= ^ 


cv'p'-p.2 


und die Laplace’sche Gleichung (/) = 0 wird 


8 /.8/\ p2-v2 8 / 8/ 


5C 8p \ 8p 


+ 


B~ 

CA 9 (ji 8 {X 




p2-|.2 0 


AB 0v 


C-^1 = 0. 

8 V 


Diese Gleichung wird aber bekanntlich durch die sogenannten Lame’schen Produkte befriedigt, das 
heißt durch 


/= RMN 
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Dr, K. Hillebrand, 


wo R nur von p, M nur von (j. und N nur von v abhängt und folgende Gleichungen befriedigt werden: 


^ ^ = [n(n+ l) R 


B — i— -\n (n + 1) p?-K^] M 
\ 2[J./ 

C-L/c—)= [n{n+\)'^^-K‘^]N. 

0 v \ öv y 

n ist eine ganze, positive Zahl und K kann so bestimmt werden, daß die betreffenden Größen ganze 
Funktionen nten Grades ihrer Argumente sind. 

MN = A ist eine Lame’sche Flächenfunktion, die offenbar der Gleichung genügt 

(p 2 _v 2 ) B -^[b^\ + fc-^] = [n (n + 1) 

OIK \ 0 [K J öv\övy 

Für das abgeplattete Rotationsellipsoid wird 


und 


daher 


b = v = 0 


lim — =: sin 9 , 
h 


dv 1 , , dii d%‘ 

hm - ~ — d^ und —^ = - 

C c B c sin ^ 


und die obige Gleichung in A, wenn sin^s gesetzt wird, 




ÖA\ 


32 A 


sin d' - sin f)- *- -h (1 — sin^ s sin^ --+ A sin^ ^ [n {n + 1) —in^ sin^ s] = 0 


0^ 1 


Bezeichnet man den Operator 

1 0 


sin O’ 0^ 


sin d' 


0^ 


0(p2v 


1 — sin2 s sin2'8- 02 

H-= V 

sin 2 0 92 


so ist 


V (A) =: ~-[n {n + sin^ s] A. 

Man sieht sofort, daß der letzten Gleichung genügt wird durch 


A = 

wo P eine zugeordnete Kugelfunktion ist. Man erhält nach Substitution für P die Gleichung 

0 / 0P\ 

sin d‘ - sin ^{n + 1) sin^^* —P=0, 

* 08’ \ 08 y 

das ist die Differentialgleichung der zugeordneten Funktion P[,f(cos8). Es ist also bis auf konstante 
Faktoren 

M = PW (cos &) = 

N = 


Da R dieselbe Funktion von p ist wie M von [a, also 
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so ist die räumliche Lame’sche Funktion bis auf einen konstanten Faktor gleich 


(i cot s) (cos) '0’ 


Gewöhnlich wird der an sich willkürliche Faktor der Lame’schen Flächenfunktion A so bestimmt, 
daß in der bekannten Integralrelation 


■ ■ ... = Konst. 

J(En.)\/(p®-!A^)(p^-v2) 


die rechts stehende Konstante der Einheit gleich wird. 

Soll nun unter die zugeordnete Kugelfunktion bezeichnet werden, wie sie in Heine’s 

»Theorie der Kugelfunktionen« definiert wird: 


pw w = 


1.3.(2«-l)^' ■ dx"^ 


{nm )! 


1)2 -i^ 

dx«+’'^ 


SO erhält man 


__ 

A izi —— \J^n + 1 

V/27: 


1.3...(2^i-l) 


\J{)i + «»)!\y’(K — in)\ 


COS e^‘ 


im 'f 


Ist Pn,mi^) Neumann definierte zugeordnete Kugelfunktion, das heißt 


Ph.H! (x) = (l-x'^) 


^ d' 


dx"' 


1.3...(2n-l) 
{n — ni )! 




so ist 


= ^ (cos &) 

\J 2iz {n + in)\ 


In der Folge wird es notwendig sein, für die von p abhängigen Teile auch die Lame’schen, 
beziehungsweise Kugelfunktionen zweiter Art einzuführen; sie genügen denselben Differentialgleichungen 
und stehen mit den entsprechenden Funktionen der ersten Art in der Beziehung 


S={2n +l)i? r 
Ja ^ 


dp 


V/ (p2-^72)(p2__^2).^2 


wo S die Funktion zweiter Art ist, die der Funktion erster Art R vom Grade n entspricht. 
Zwischen beiden besteht die Relation 


dS d R 2 n \ 

'di~ 


Die entsprechenden Relationen für die zugeordneten Kugelfunktionen (nach der Heine’schen 

Definition) lauten 



und 


n«) (.r) 


dQ^y) 

dx 




dP^^ix) 

dx 


2n + \ 
1-x^ 
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Nach der Neumann’schen Bezeichnungsweise ist 

ö„, (^) 2 ^)! gw 

und die entsprechenden Relationen 

Ö„, (X) = 2 p„^ (;,) f 


dx 


Pn,,n(x)-^^f^ - Qu,,nix) ^ 

dx dx * 


{l-x^)iP„,,„ix)Y 
dP„,m(x) _ 2 (n+m)\ 1 


(n — m )! l~x^ 


Arten der zugeordneten Funktionen 


Um für das abgeplattete Rotationsellipsoid die Funktionen R und S zu finden^ hat man in den beiden 

X ~ . i — = i cot e 

V 


zu setzen, so daß 


und die Differentialbeziehung 


wird. 


dx 


cd^ 


1 — x^ p\/c^—p^^ 


z= —id 


p(„) ^ _ g(„) ^ (2w + 1) 

ÖS as 


4. Die beiden Hauptgleichungen der Theorie der Gezeiten. 


Vermöge der Darstellung der beiden Unbekannten n und v durch die Funktion ^ werden die 
Bewegungsgleichungen identisch befriedigt und es handelt sich nur darum, diese Funktion ^ und die noch 
weiter dazutretende unbekannte Störung C der Niveaufläche so zu bestimmen, daß der Kontinuitäts¬ 
gleichung und der Oberflächenbedingung genügt wird. Vorausgesetzt ist in der oben angegebenen Form 
schon ein ganz bestimmter funktioneller Zusammenhang mit der Zeit, demgemäß die Verschiebungs¬ 
komponenten sowohl als auch ^ und C nach Exponentialfunktionen entwickelt gedacht sind, die im 
Falle der erzwungenen Schwingungen den Gliedern in der Entwicklung der Störungsfunktion entsprechen, 
während im P'alle der freien Schwingungen sich aus den Bedingungen des Systems die X bestimmen lassen. 

Da es sich bei den vorliegenden Betrachtungen nur um periodische Vorgänge handeln kann, so soll 
in der Folge /X statt X eingeführt werden. Die in der Kontinuitätsgleichung vorkommenden Größen und 
/Zg sind dann 


h = 


h 


1 _ iiü! cos« 6 
X« 


^2 =- cos 

Z X 


Da ferner für das Rotationsellipsoid 

E p^ (1— sin^ £ sin^ '9’) , G = p^ sin^ 9’ 


ist, so hat man als Kontinuitätsgleichung 



_ ^ _ ^ fj^ \/1 — sin^ s sin^ 9^ ^ ^ \ 

—sin^ s sin^ 9* d d’J 8cp\ sin9' Bcpy 

2(0 /d {h^ cos 9*) 8 (j; 8 cos 9-) 

A \ 89- 8 cp 8 cp 



= sin 9’\/l — sin^ sTin^ 9^. C 
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und als Bedingung der freien Oberfläche 

-^0 v/l—sin^ssin^a'.C + +5=0. 

Die Lösung der Aufgabe wird nun dadurch bewerkstelligt, daß man sich die Unbekannten ^ und 
\/l—sin^ssin^^.C als Funktionen des Ortes auf der Gleichgewichtsoberfiäche entwickelt denkt nach den 
dem Rotationsellipsoid zugehörigen Flächenfunktionen P^lf^(cos '9-) und aus den beiden Bedingungs¬ 
gleichungen die Koeffizienten dieser Entwicklung zu ermitteln sucht. Dann sind die Ableitungen unmittel¬ 
bar aus der Natur dieser Entwicklungsfunktionen gegeben und die Lösung des Problems wird sich auf die 
Auflösung einer Serie endlicher Gleichungen in den Koeffizienten reduzieren. 

Zunächst sieht man, daß für jedes dieser Glieder 


—^ ~ 
d (p 


02 ^ 

—I. 6 

8 cp 2 


ist. 

Es soll nun weiter die Annahme gemacht werden, daß die Tiefe der ungestörten Wasserschichte 
eine zonale Funktion ist, so wie es bei völlig exakten Niveauschichten der Fall wäre, und Abweichungen 
davon einer späteren Untersuchung Vorbehalten bleiben; dann ist h^ nur von d' abhängig und die Kon¬ 
tinuitätsgleichung nimmt die Form an 


/ K sin ^ _ /^yi-sin^ssin^^ ^ ^ ^ 

8'0' \\/l — sin2 s sin2'0' sin 0 X 8 0’ 

= p2 sin 0’'^!— sin2 s sin^ 0-.^. 

Dieser Gleichung kann man eine etwas andere Gestalt geben, die für gewisse einfache Annahmen 
eine besonders leichte Behandlung gestattet. 

Setzt man der Kürze halber 


so daß 


2 CO 

T" 


= 0, 


h 

1 —a2 cos2 0 


führt statt dieser Größen h ein 


H = - 


\/l—sin2 s sin2 0- 
und bezeichnet mit dem Operator 


, H,= 


K 


\/1 —sin2 s sin2 0* 


D = — sin 0* 


80- 


so erhält man nach einigen leichten Transformationen für die Kontinuitätsbedingung schließlich die folgende 
Gleichung 

(D + om cos 0') [H^ (D — om cos 0’) (]>] ~H (1 —sin^ s sin2 0’) tj) 

= p2 sin2 0'N^ 1 —sin2 s sin^ 0*,C’ 

Was ferner die Bedingung der freien Oberfläche anbelangt, so kann zunächst das darin vorkommende 
Integral 
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durch die Entvvicklungskoeffizienten von C \/l — sin^s sin ^0- selbst erhalten werden: es ist ja nichts 
anderes als das Potential einer unendlich dünnen Schicht C auf einen ihrer Punkte und als solches eine 
bestimmte Funktion des Ortes auf der Oberfläche des Rotationsellipsoides, kann also wieder nach Kugel- 
flächenfiinktion d' und cp entwickelt gedacht werden, und zwar soll das in der Form geschehen 

n = SaW/p(«) _ 

wo für den Moment der Äquatorhalbmesser der ellipsoidischen Schicht gesetzt wird, zum Unterschied 
von einem veränderlichen p. 

Dann ist das Potential dieser Schichte auf einen außerhalb gelegenen Punkt mit den elliptischen 
Koordinaten p > po, cp 

n,= SaW 1 - PÜlHcos ») .'-f 

und das Potential auf einen inneren Punkt p <Pq, cp 

n,- zu Sa(«)/P«) Pi«>(cos , 


A. 


0(’‘) 


1 - Po 


PW (cos ■ö-) ef 


m 'f 


denn die Reihen konvergieren zugleich mit 0, genügen der Laplace’schen Gleichung = 0 und nehmen 
auf der Oberfläche selbst die vorgeschriebenen Werte an, stellen also nach dem Dirichlet’schen Prinzip 
tatsächlich das äußere und innere Potential der Schicht dar. 

Die Koeffizienten a sind offenbar reelle Größen. Nach der Heine’schen Definition der zugeordneten 
Funktionen ist ' • 

PW (i cot s) = 4 (i cot s) = SW, 

wo 


R(n) _ 

cos”'“'” S 

P 

(n “ m) (n — m — 1) 

sin ”s 

2.(2«-l) 

= 

— 

sin ” + U 

h- 

(n + m + 1)(^^ + m+2) 

cos”+'”+^ s 

i 

2.(2«+3) 


(n — nt). . . (n ~ m — 3) 

2.4.(27^-1) (2^;-3) 

{n + 77^ q- 1). . , (n + m + 4) ^ 

2.4.(277 + 3) (277+5) ^ ^ 




■] 


Es besteht dann nach der obigen Relation für die beiden Arten der zugeordneten Funktionen die 
Gleichung 


^ -Sä) - 2n + 1 


8s 


8 s 


und 


n - S a(«) (cos 

(cos^)^^-^»'^ 

n, = (cos^)^^'-'^, 

wo der zweite untere Index 0 bedeutet, daß das Argument für p^ zu nehmen ist. 

Ist weiter die Dichte der Flüssigkeitsschicht und betrachtet man C konsequent als kleine Größe, 
so kann die zu IT gehörige Flüssigkeitsmenge als unendlich dünne Schicht mit der Dichte Sq C angesehen 
werden und dann ist nach einem bekannten Theorem 
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(d^a dUA . ^ 

lim - - -^1= — 4x5o C, 

\ du dnJ 

wo die Ableitungen nach der nach aussen gerichteten Flächennormale zu nehmen sind. Aus dem Aus¬ 
druck für das Linienelement in elliptischen Koordinaten folgt aber 


demnach 


dn — 


1 

cos e 


\/l — sin^ s sin^ -d- dp 


und wegen 


4xSoC 


— 4 7C 5o C 


cos e f dUa 

, ■ —- - - hm -—- 

— sin^ £o sin^ O’ \ dp 


äUj \ 

^ P J9-9o 


—A — — cot s de 
9 


1 sin So 

pQ s/l— sin^ So sin^ 'ö- 





d 

ÖS 


ö \ 

^ P(”)(cos^) 

ÖS j 


und infolge der obigen Relation, wenn wieder p für p^ gesetzt wird, 


5o C v/1 — sin^ s sin^ & — 


^ E(2n -h 1) Pj^^cos d’)e^'^'^'^ 

47rp 


woraus hervorgeht, daß sich die Koeffizienten der Entwicklung des Potentiales 11 unmittelbar durch die 
der Reihe für 

C\/l — sin^ s sin^ & 

ausdrücken lassen. Ist die Entwicklung für letztere Größe gegeben durch 

C \/l — sin^ s sin^ 0- =: S^^5)p(")(cos 0-) 
so erhält man für dieses Potential 

n = P^o . 

^ 2n -h 1 sin s 

Die Bedingung der freien Oberfläche ist demnach 

sin^s sm-^^.C-X2(];-Fy^2n + S = 0 


und gibt zusammen mit der Kontinuitätsgleichung die Lösung des Problems, die also darin bestehen wird, 
die Koeffizienten der Entwicklung der beiden Unbekannten C und ^ nach Kugelflächenfunktionen von ^ 
und cp zu bestimmen. 

Denkt man sich also — immer für ein bestimmtes X und abgesehen von dem Faktor e^'^^ — 


daher 

wo 


C s/1 — sin^ s sin^ O* = S Pjjj) (cos d') 


n = A(f) P(«) (cos , 


Jr7n — 


4ir 


2m + 1 


sin s 


also eine bekannte Funktion von s ist, ferner 
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D\\ K. Hill eh r an d, 

fh = (cos , 

so wird es sich um die Ermittlung der Größen A^H^und T^JJ) handeln. Setzt man den zu X gehörigen Teil der 
Störungsfunktion 5 ebenso entwickelt voraus, also 

S = S Q/) (cos , 

wo die bekannte Größen, im Fall der freien Schwingung gleich Null sind, so müssen nach einer 
bekannten Eigenschaft dieser Entvvicklungsfunktionen nach Substitution in die beiden Gleichungen die 
Faktoren derselben Kugelfunktion indentisch verschwinden, eine Bedingung, die eine Serie von Gleichungen 
zur Bestimmung der Koeffizienten und ergibt. 


5. Methode der Lösung für den Fall einer zonalen Tiefenfunktion. 

Nach dem eben angegebenen Lösungsvorgang erhält man aus der Bedingung der freien Oberfläche 
die die eine Serie von Bestimmungsgleichungen repräsentierende Relation 

Weniger einfach gestalten sich die aus der Kontinuitätsgleichung folgenden Bedingungsgleichungen. 

Das Nullsetzen der Koeffizienten der linken Seite der Gleichung wird dadurch bedingt, daß die totale 
Funktion, die diese bildet, nach Kugelfunktionen entwickelt ist. Nun kommen hier die Größen und 
C\/1 — sin^ 8 sin^ d', für welche diese Entwicklung supponiert wurde, in Verbindung mit gewissen Funktionen 
von “ö- und mit dem Operator D vor. 

Es sind also zunächst diese Verbindungen der in und C\/l—sin^ s sin^'9’auftretenden Kugel¬ 
funktionen wieder durch solche Funktionen auszudrücken. Das ist tatsächlich möglich. Aus der bekannten 
Formel für die zugeordneten Funktionen 


^ ^ + nP(f/+^Xx) ~ 
dx 


(n — m) {n H- m) 
{2n-\){2n+\) 


(h +1) = 0 


folgt unmittelbar, daß 


D (PW (cosd)) = - Mi^»+i)(cosd) + (« + 1) (cos d). 

(2»—1) (2n + 1) 


(Die Heine’sche Definition vorausgesetzt.) 
Ferner folgt aus der Rekursionsformel 


;rPi«)(^) = P<^+i)(A-)+ 


(n+m) 

(2n-h l)(2n-l) 


und deren wiederholter Anwendung, daß jedes Produkt 

cos^'^.P(JJ)(cos^) 

sich ausdrücken läßt durch 


p(n~\-k) p(n~\-Ji—2) p(n~k) 

wenn k ^ n, respektive durch 


wenn li> n ist. 


P(«+^) p(w+^~ 2 ) p( 0 ) 

^ m y m y” ' ' ” ^ my 
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Denkt man sich das Tiefengesetz in einer Potenzreihe nach cos i> hergestellt, so wird man die 
Kontinuitätsgleichung in Kugelfunktionen mit konstanten Koeffizienten auflösen und so die Bedingungs¬ 
gleichungen für diese erhalten können. Es werden dies Relationen zwischen einer endlichen Anzahl von 
Koeffizienten und FJ”) sein, mit gleichem unteren Index m und gewissen Reihen aufeinanderfolgender 
oberer Indizes. 

Es hat nun keinen Zweck, diesen allgemeinen P"ali weiter zu verfolgen, der ja auf ganz willkürliche 
Konfigurationen Bezug nehmen müßte; man wird, um allgemein gültige Näherungen zu erhalten, einen 
gewissen mittleren Verlauf der Tiefe nach der geographischen Breite annehmen, für welchen die Lösung 
sich wesentlich einfacher gestalten wird. 

Es sei hier nur noch ein Laplace'sches Theorem erwähnt, das sich, wie Poincare bei dem ein¬ 
facheren Fall der sphärischen Gleichgewichtsoberfläche bemerkt hat, unmittelbar aus dem Anblick dieser 
Lösungsform ergibt. 

Aus dem ganzen Vorgang der Lösung ist ersichtlich, daß, wenn es sich um erzwungene Schwingungen 
handelt, jedem Glied der Störungsfunktion, das e^'^^ als Faktor enthält, ein Glied mit dem gleichen 
Faktor im Ausdrucke für und C\/l—sin^ £ sin^ ^ entspricht. Der übrige Faktor ist eine gewisse durch 
die {nt Konst.) ausgedrückte Funktion von d', deren Koeffizienten sich eben aus dem Gleichungssystem 
bestimmen, die aus den beiden Hauptgleichungen erhalten wurden und daher sc wie diese nur reelle 
Größen enthalten; es werden also auch die zu ermittelnden Koeffizienten entweder reelle oder paarweise 
konjugierte Größen sein und — wenigstens in ihrer Vereinigung — keinen Faktor von der Form e^^ liefern, 
wo k irgendeine Konstante bedeutet; das heißt aber nichts anderes als daß die entsprechenden Glieder 
der Gezeitenbewegung dieselbe Phase besitzen wie die störenden Kräfte. Das ist der Inhalt des einen 
Laplace’schen Theorems, welches besagt, daß keine Flutverzögerung eintritt, wenn die Tiefe der Flüssig¬ 
keitsschicht nur eine Funktion der Breite ist. 


6 . Lösung für eine exakte Niveauschichte. ^ 

Setzt man — was bisher immer stillschweigend getan wurde — eine vollständige Bedeckung des 
festen Himmelskörpers mit der Flüssigkeitsschichte voraus, so wird man gewisse mittlere Verhältnisse 
erhalten, wenn man die Unregelmäßigkeiten der festen Oberfläche vernachlässigt, demnach ein solches 
Tiefengesetz annimmt, wie es stattfinden muß, wenn die Flüssigkeitsschichte zwischen zwei unendlich 
benachbarten Niveauflächen eingeschlossen ist: Aus der Beziehung zwischen dem Element der Flächen¬ 
normale und dem Inkrement des Potentialwertes der Schwerkraft auf der Oberfläche folgt, daß die Tiefe 
umgekehrt proportional der Beschleunigung der Schwere sein muß, also 


1 — sin^ £ sin^ ^ 

wo h^ eine Konstante ist, eine Annahme die der konstanten Tiefe bei sphärischer Gleichgevvichtsober- 
fläche äquivalent ist. Es nimmt dann die Kontinuitätsbedingung eine wesentlich einfachere Gestalt an. 
Sie wird am leichtesten durch eine Transformation erhalten, die eine Verallgemeinerung einer von 
S. S. Hough in der eingangs zitierten Abhandlung »On the Application of Harmonie Analysis to the 
Dynamical Theory of Tides«, Part. II, Phil. Trans. A., vol. 191, angegebenen Umformung dieser Gleichung 
für eine sphärische Oberfläche ist. 

Substituiert man für zwei neue unbekannte Funktionen und vermöge der zwei Relationen 
(j> = (Z) -P am cos^) -P (1—sin‘^ £ sin^^) (1—a- cos^6)({j.^, 
sin^ö* (V -P am) Z) [(1 —sin^ £ sin^ ö-) (1 — a' cos- 6)], 
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Dr. K. Hill ehr and, 


wo V der in Nummer 3 definierte Operator ist und die zweite Gleichung mit Berücksichtigung der Be¬ 
deutung von D auch 

(V+ om) — ~ 2 cos ^ (sin‘^ s —o^) 

geschrieben werden kann, so nimmt die Kontinuitätsgleichung die Form an 
(D + 0 m cos {H (1 — sin^ s sin^ [ni^ + (D — g m cos 

— H.( \ “Sin^ e sin^ [(D -h g m cos i>) 4- (1 — sin^ s sin^ '8-) (1 —a^ cos^ 6) 


= C sin^ ^ 1 — sin‘^ s sin‘^ d '. 

Macht man die Annahme, daß 

H(1 —sin^ s sin^ ^)r= 

einer Konstanten gleich ist, so fällt aus der Gleichung heraus und die Kontinuitätsbedingung reduziert 
sich aut 

K (V —am) = Cs/l^-sin^ s sin‘^^. 

Die obige Annahme ist aber nichts anderes als 

h 

h = 


\/1 — sin^ s sin^ ^ ^ 

das heißt, jener Verlauf der Tiefe, der stattfindet, wenn auch die feste Oberfläche eine Niveaufläche ist. 

Zur Lösung des Problems hat man jetzt vier Gleichungen in den vier unbekannten Funktionen 
C, und 

Berücksichtigt man in der Darstellung von ^ durch und daß 

cos-ö- 


cos 6 — 


so lauten diese 


\/1 — sin^ s sin^ 


— gQ \/1 —sin^ £ sin^^. C —^ + S “ 0 


7^0 (V—am) = p^C s/1 —sin^ s sin^ d- 
(V + am) = — 2 cos ^ (sin^ s*—a^) 
r= (Z) 4- am cos ^) + [cos^ s + (sin^ s —a^) cos 2 ^](j^ 2 . 


oder 


(J> = (Z) + a m cos 'S-) -cos -ö- (V + am) + cos^ £ • ^2 • 

2 

Die Lösung geschieht nun so, daß man auch für die Hilfsfunktionen und Entwicklung 

nach Kugelflächenfunktionen voraussetzt, im ganzen also vier Reihen von Koeffizienten aus den vier 
jetzt wesentlich einfacheren Gleichungen zu bestimmen hat, und zwar sollen hier die Neumann’schen 
Zugeordneten eingeführt werden, die eine bessere Symmetrie in den einzelnen Relationen ergeben. Für 
diese lautet die hier benützte Rekursionsformel 


und 


{2n 1) cos ^ Pn^,n (cos f)*) =(n — m + 1) Pn+I,m (cos &) + {n + m) Pn-x,m (cos ^) 

{2n 4 - 1) D[Pn,,n (cos ^)] — — n {n-~m -t- 1) Pn-\-i,m (cos ^) + (;^ — 1) {n + m) Pn-i,m (cos -O-). 
Es sei also der Faktor von 
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in 

C\/l — sin^ s sin^ ^ : hAn^m Pn,m (cos tl) 

daher in 

n : An, nt Pn, in(cOS 

Pn,ni (cOS 

^2 • ^?n,m Pn,m (cOS -ö*) 

S . ^Cn,m Pn,m ' 

Führt man diese Reihenentwicklungen in die vier Gleichungen ein, berücksichtigt außer den obigen 
Rekursionsformeln noch^ daß 

V [P«,m(cos^)] = — [n(n + l)—m^ sin^ s] Pn,m (cos 9-) 

und setzt die Koeffizienten der einzelnen Kugelfunktionen gleich Null, so erhält man die folgenden 
Rekursionsformeln in den zu ermittelnden Größen A, F, a, ß, wobei der Einfachheit halber der allen 
gemeinsame Index ni vorläufig nicht eigens angemerkt werden soll: 

A, + \^Vn=Cn 

h^l—n {n + 1) + f/P sin- e—am] ß„ zu An{A 

[— n (n + 1) + sin^ s -f- oni] = 2 (a^ —sin^ s) ß„_i + nt -h 1 \ 

\ 2^/—1 2;/ + 3 / 

1 

r„ cos^ £. + — a«-i- [0^— 1) (n~2) — nf-‘ sin^ s 4- a ni\ 

2 27^-1 

4 —— y.n +1 - - 10^ + 2) 0^ + 3) —77^^ sin^ s + am]. 

2 277 + 3 


Durch Elimination der a, ß und F erhält man für die A die folgende Rekursionsformel 

^(32_sin2e) («-^^-1) («-^^) _ ^ _ 

h^ (2n — S)(2n—l) (n—l)n — atn—m^ sin^ s 




^ _ 

Jiq 77 ( 77 + 1 )- 


cos^ s 


-am 


m^ sin'-' s 


+ 


+ 




a*^ — sin^ e 


1) —0 777 + sin^ s 


(n — m) (77 + 777 ) (77 — 2) (77 — 1) + am — sin^ e 


(2 77 — 1) (277 + 1) (^^—1) om — m'^ sin^ s 


(77— m + 1) (77 + 777 + 1) (77 + 2) (77 + 3) + am — m^ sin^ ^' \ + 

(277 + 1) (277 + 3) (77 + 1)(77 + 2) —0777 —7772 sin‘^ s _ I 

/ 2 sin^ s) (^^ + 777 + 1) (77 + 777 + 2 ) _ An -[.2 _ _f+ Cn 

JtQ (277 + 3) (277 + 5) (77 +1) (n + 2) —0777 —777^ sin^ £ ( 0 


Führt man folgende Bezeichnungen ein 

a, = ^ („>_ sin- s) fc"‘ + 

h 


1 


(277+ 1) (277 +3) (77 + 1) (77 + 2) — 0 777 —777^ Sin^ £ 
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/ o • o ^ (^ + « 2 — 1 ) (n + m) 1 

e„ = —^ (a^ —sin® e) - -^--- 

kg (2 22— 1) (2 22 + 1) (22 — 1) 22 ^ 0 222 — 222® sin® S 


wo 


und 


Bn = 


1 


x®p® _ 

hg n (n + 1) + 0 222 


— Sn + B„, 

in — 

- / — cos^ s + 

— 7;^2sin2£ ( 


(a- —sin^ s) 


{n — 'iu) {n + m) {n — 1) {n — 2) + om — m- sin^ s 


+ 


{2ti— 1) (2«+ 1) {it— 1) n — m^ — ni^ sin^ s 
1) {n H- ni + 1) {n +2) {n H- 3) + ain — m^ sin^ s 


+ 


(2n + l){2n + 3) (// + 1) (n + 2) — am — m^ sin^ s 

so lautet die Rekursionsformel 


an~2 An-2 + 


= 

\o ’ 


Dabei ist aber zu bemerken, daß wegen 

Pn,m (cos^) = 0 für n<m 

diese Formel überhaupt nur für n ^ m gilt und in den ersten beiden Fällen n ~ in und n — in -p 1 

■^m—2 — -^m—l — ^ 

ZU setzen ist. 

Aus der Bedeutung von C folgt ferner wegen der Unveränderlichkeit der Quantität der deformierter 
Masse, daß 

J(:do = 0 

sein muß, wenn das Integral über die ganze Oberfläche ausgedehnt wird (immer vollständige Bedeckung 
mit der Flüssigkeitsschicht vorausgesetzt), das heißt also, 


oder 


J J Csin'G’Syi —sin^ s sin^'ö^.^^fl-t^rp — o 


JJSA.h! PH,,„(coS'&)e”''»sin = 0. 


Da vermöge der bekannten Eigenschaften der Kugelfunktionen diese Gleichung durch jedes Glied 
der Summe für sich erfüllt ist, außer für m =: n = 0, so muß 


^03 0 — ^ 


sein. 


Man kann übrigens dem Koeffizienten Bn eine etwas andere Form geben, die die Behandlungsweise 
der Bedingungsgleichungen erleichtert. 

Zerlegt man-die von aabhängigeir Teile in Partialbrüche, setzt man - - 

(n — 1) (n — 2) —m^ sin^ s + cm 


[n (it + — sin^s + cm] [{n—l) n — m^ sin^ s — om] 


1 r 2n-l 


- - f. I C 


{n~ 1)2 —sin2 s 


n {fl + — smH + am {n—\)n — m-s\ri^z—om 
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und 


{n + 2) {n + 3) —sin^ s + om 


[n {n +1)—s+a;;/] [( 7 / + 1) (77 + 2) — sin^ 


1 


( 7 ^ + Vf — m^ sin^ s 


277 + 3 


(77 + 2)^ — 7772 sin2 £ 


77 (77 + 1) — 777^ sin^s + (77 + 1) (77 + 2) —777^ sin^ z — 


so erhält man 


X^p^COS^S 77 (77 + 1) —(3 777 --7772 sin‘^ £ 

Dn — — -- 


[( 77 + 1)2 — 7772 sin®s] [ 772—7772 sin2 s] 


X2p2. 


0^ —sin^ £ 


{11 ~m) {11 + 777) [(77-- 1)2 — 7772 sin2 s] 


(277— 1) (277+1)(772 — 7772 sin2 £) [(77— 1) n — om~ni^ sin2s] 


(11 —m + 1) (77 + 777 + 1) ([77 + 2) 2 — 7772 sin2 s] 


(277 + 1) (277 + 3) [(77 + 1)2-7772 sin2£] [(77 + 1) {n + 2) 


n^s] , ) 

-~(sm~nT^ sin2 b]\ 


Es soll hier noch eine Umformung angegeben werden, die sich besonders bei der Ermittlung der 
erzwungenen Schwingungen, also bei gegebenen X als vorteilhaft erweist. Die Gleichung 

a n- 2 Än —2 + ^«+2 ~ Cn 

kann nach Substitution der angegebenem Werte für die + h und c auch geschrieben werden 

77 ~ 7/7 


r 


(277— 1) [(77 — 1) 77 — 0777— 7/72 sin2 £] ( 2/7 — 3 

Jh 

X2p2 


— 7/7 — 1 ^ ^ (77+7/7) [(77— 1)2 —7/72 sin 2 s] ^ 

^«-2+-^^-1—;— 


+ 


■ sm2 £ ( 


(2 77 + 1) [772 — 7/72 sin2 £] 
2 ein 2 5 . 


L -C0S> + - 

[( 77 + 1)2 — 7772 sin2 £] [772 — 7/72 sin2 £] 

77 + 777+ 1 f/7 + 7/7 + 


1 


(277 + 3) [(77+ 1) (77+2) —0 777 —7772 sin2 £] ( 277 + 5 

//q Cji 

X 2 p 2 o 2 —sin 2 £ 

Führt man nun Hilfsgrößen D ein vermöge der Relation 

(77 —777) [(77+1)2 —7772 sin 2£] (77+7/7+ 1 ) (7/2 —7/72 sin2 £) _ 

Aji-i + -::-- ^«+1 — 


+- 2 . (77 — 7/7 + 1 )[(77 + 2)2 — 7/72 sin2 £] ] 

- A}i^2 H---- An} 

5 (277+1) [(77+1)2 — 7/72 sin2£] j 


2/7-1 


so wird aus der obigen Gleichung 


277 + 3 


= [77 (77+l)—o7/7 — 7772 sin2£] Dn, 


n—m 

277— 1 


[(77+ 1)2 — 7772 sin2 £] Z)«-!" 


gnh 


— sin2£) [(77+ 1)2 — 7/72 sin2 £] 


— C0S2£ [77 (77+ 1) — 0 777 —7772 sin2£] 


An +1 9•2^n 

-j- ^— sin2 £) Dn+x = 


o2 —sin2 £ 27^ + 3 


_ h ^ — sin2 z ) [(77 + 1)2 — 7/72 sin 2 s] 

— . . 

X2p2 ^ ^ 


0^ —sin^ £ 


Setzt man 

77— 7/7 + 1 


[(7/+ 2)2 —7/72 sin2£] = ktiy 


[(77— 1)2 — 7/72 sin2 £] = 

2/7 + 1 


2/7+1 
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sin^ s 


Icos^ e [n {n-^ \ ) — om — iit^ sin^ e]- 


K gn 

X2p2 


(n^ — m- sin^ e) [{n-^iy — m^ sin^ b]\ = 


so ist der Gleichung 
äquivalent das System 


Jh_ 

X2p2 


n (n-hi) — oin — m- sin^s = Nn 

sin2 s) [(;/-!- 1)2 _^;^2 sin2 s] ^ 

—sin^ e 

an-2 + + ^;j-f2 = Cn 

^«—1 Ai—1 — ^« + 4+1 Ai- fl = Kn 

kn—2 II—2 71—1 Ai—1 ”1“ hi ~ ^ 


und die beiden Systeme 

^iji -^7» ”b -^7«+2 ~ Cni 

Cljn -^771 ~l~ ^;7i-f2 -^m-f2"j”G;i-f4 -^m-f4 — Cjn~\-2 


— Ai- fl -D^-f i + ln+2 ^fi-f 2 — 0 


^771+1 -^7;i-fl "1” ^711-1-3 -^771-f 3 - Al-fl 

^771-fl -^»I-f 1 H" ^771-f8 '^77I-f3“b ^77lf-5 -^77I-f5 - ^77I-f3 


werden ersetzt durch 

^771— i£> 771—1 -^771 -^771 4" A-f 1^771-fl K^l 

Al-flAll-fl —A71-f2^771-f2 +A-f3 ^771+3 = A71-f2 


Al ^771 -^^Ai-f 1 4“ lm-\-2 -^771+2 — -Aii-f 1 

^771-f2-^771-f2 ^Al-f3 -^771-f 3“l“^771-f4 ^771-f 4 ~ A71-f3 


A—1 Dni —1 + Int ^771 — ö 
^771 -^771 ^7n~\-l ^771-f 1 “l“ A-f2 -^771-f2 — ^ 


A?;^ Dfn~^ ^771-f 1 -^771-fl - 0 

^771-fl^771-fl—A;„^2A71-f2H“A-f3^771-f3 = 0 

! 

) 


7. Freie Schwingungen. ^ 

Da für die freien Schwingungen sämtliche C = 0 sind, so sind die Koeffizienten nach dem j 
bisherigen bestimmt durch zwei Systeme rekurrierender — im Allgemeinen trinomischer — Gleichungen | 
von folgender Form, wobei wieder der allen Größen gemeinsame zweite Index in vorläufig unterdrückt 
wird: 

^7M ^771 + G71-f 2 ^77I-f 2 ~ 0 ^ 

^771 -^771 ”F ^771-f 2 -^»l-f 2 “F G7I-f4 -^771+4 0 

^77I-f2 -^771-f2 “F Al-f4 -^771-f4 “F Cin-\-6 -4772-f6 — 


i 


a,j_4 An—4’^hn~2 ^7i—2 + Gi ^ii — 0 
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4“^w-|- 3 4“ -^m+ö ~ ^ 

^w+3 ^w+3 + ^w+5 ^m+5 + ^m+7 ^w+7 — 0 


^«—3 An—Z + bn~\ An—\ + CnJ^\ AnJ^\ — 0 . 


Das erste System gibt die Koeffizienten jener Pn,m (cos '8-) für welche n — m eine gerade Zahl ist, 
das zweite für solche mit ungerader Differenz n — m\ das erste definiert der Bedeutung der zugeordneten 
Kugelfunktionen gemäß Schwingungen, die symmetrisch zum Rotationsäquator stattfinden. 


Die Lösungen derartiger Systeme können bekanntlich in Form von Kettenbrüchen dargestellt 
werden. 

Setzt man 


An-~2 j 

^n-~2 —:— — nn~2 , 

A„ 


^n-^2 



1 


SO findet man aus den beiden Gleichungen 


und 

daß 


^«-2 4-^«+ ^«-1-2 — 0 

Un h^i ^w-f 2 9 


an Cn-\.2 

Ky2 = — - - 


^?i+2- 


^«+2 ^w+4 


^«+2v 4- ^n+2v+2 


hn 


Un Cny2 


hfl - 


an~2 

hn-2~ 


2 v 4 -/ 2«—2 v _2 

der zweite Kettenbruch ist ein endlicher, weil für n < m 


An, m — 0 


ist. Hat man auf diese Weise die Größen h und^ ermittelt, so werden sich sämtliche An der ersten Gruppe 
durch Anti die der zweiten durch Ain-\-i dargestellt ergeben. 


Damit aber diese Größen von Null verschiedene Werte erhalten, tritt hier noch die Bedingung dazu, 
daß die Determinante jedes Systems verschwindet. Diese kann für das erste System in folgender Weise 


formuliert werden — für das zweite gelten analoge Beziehungen —: setzt man die Determinante 
Grades 


n—m 


2 


1 Eine Verwechslung mit den k der vorigen Nummer ist ja nicht zu befürchten. 
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0 

0 

0 


^ 0 0 

^;«+2 0^ 0 

^ ^m+6 0 


^«—6 —4 0 

Uji —ß ^«—4 Cji —2 

0 a^i^^ bfi—-2 


so muß 


lim Dfi — 0 


n—oo 


sein und das ist offenbar eine Bedingung für die in den a, b, c enthaltene Größe \ demnach die Gleichung 
aus der die Perioden der freien Schwingungen der den Himmelskörper bedeckenden Flüssigkeitsschichte 
erhalten werden. 

Die Lösungen X werden eine nach einem bestimmten Gesetze fortschreitende unendliche Reihe 
bilden, jeder Lösung gehört ein bestimmtes Wertesystem 



an, worin also A^u beziehungsweise öffenbar die willkürlichen Konstanten des Problems sind. 

Für die numerische Auswertung der X und der zugehörigen Koeffizientensysteme A^ : Aut wird sich 
nun ein Weg durch folgende Überlegung finden lassen. 


n—in 

Zur Ermittlung irgendeines A^ aus der willkürlichen Konstanten Aui genügt das System der —-— 
Gleichungen 


Aui-\-2 

^m-\-2 Aui-\-2~^l^m-\'4^ -^/«-|-4~PG«-|-6 


aui^A ^;«+4H"^W:f6 Aui-\-6-^Cui-\~S 



a„—4 An—4-^lhi~2An—2-^Cu A^ — 0. 


Da die Determinante dieses Systems 


Gr^2 


ist, so ist unmittelbar ersichtlich, daß^„ G;i-j- 2 - • 



0. -b,nA„i 

0 Ulli A}ii 

0 0 


0 


0 


—4 —2 ^ 


oder 


n~m 
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ist, woraus folgt, daß, wenn die \ der obigen Bedingung genügen, die gegen Null konvergieren. Man 
wird also jedenfalls eine Näherung für die Ermittlung der X und A erhalten, wenn man die Reihe der A 
bei irgendeinem genügend großen n abbricht, das heißt, 

An -\-2 — ^«+4 —.— 0 

setzt. 

Die Periodengleichung ist dann 

Djiy2 = 0 

und wird eine endliche Zahl von Lösungen X ergeben. 

Auch für diese läßt sich ein Näherungsverfahren angeben. Bedenkt man, daß die Größen Un und Cn 
mit wachsendem n gleichfalls gegen Null konvergieren, und vernachlässigt man schon c«, so wird 

T)n-\-2 — Gii 

und es werden in erster Näherung gewisse Werte von X offenbar aus 


erhalten. 

Man wird also aus den Gleichungen ^ 


hn = 0 


hm — 0 7 h}ji^2 — 0 y... .bji — Ö. . . . 


die Reihe der Werte für X mit einer gewissen Annäherung finden, die umso größer sein wird, je größer n 
ist. Für die Werte aus den ersteien dieser Gleichungen ist dann allerdings eine nochmalige Durchrechnung 
und daraus sich ergebende Verbesserung notwendig, das numerische Verfahren zeigt aber — wenigstens 
für den F'all unserer Erde — eine genügende Konvergenz, um uns über die Natur dieser Oszillationen in 
Kenntnis zu setzen. 

Die Wurzeln der Gleichungen = 0 lassen sich aber auf ziemlich einfachem Wege ermitteln. 

Für m — 0 sind sie ja unmittelbar anzugeben. 

Ist m eine von Null verschiedene Größe und denke man sich die X als Abszissen eines ebenen recht- 
winkeligen Koordinatensystems, dann gehören die verlangten Wurzeln den Durchschnittspunkten einer 
Kurve 

y—~Bn 

und einer zur X —Achse parallelen Geraden 
an. Die Kurvengleichung hat die Form 




Sie hat also zwei zur j'-Achse parallele Asymptoten 


A^z \ 




X= „nJX = A 


und eine parabolische Asymptote, deren Gleichung von der Form ist 


y = a'X^ — b\. 


Die Kurve wird demnach vier Arme im Raum der positiven besitzen, daher vier reelle Wurzeln 
ergeben. Für relativ große gn — deren relative Größe ja von der Annahme über h^^ abhängt - werden die 
beiden Werte 



X = 


h 
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offenbar Näherungen bedeuten. Die Wurzeln sind ganz allgemein eingeschlossen zwischen den Werten 


— oo, 0, 


2 


(die Kurve geht ja durch den Nullpunkt). 

Mit diesen aus den “0 gefundenen Näherungen wird man aus 

die verbesserten, beziehungsweise die als definitiv anzusehenden Werte der \ finden. 

Für die numerische Auswertung soll s als kleine Größe betrachtet und nur sin-s mitgenommen 
werden. Dann ist 


und 


2 . \ 

P'.. =:-I 1 H-Sin- 3 

V'- \ 5 




J 


P 


(«) 


4x ' / (2/7/-l)(2;n4-l) 

- . i + 


2/7+1 


'I 


2 (2//~l)(2n + 3) 


Sin- s 


Ist s die mittlere Dichte des ganzen Himmelskörpers, so bildet sich für 


4:: 

gn — 


<b» — <b0 

I 1 . 2 3 

1-sin^s 

10 2;/H-1 5 


'1. 


(2in—l)(2in-hl) . , 

-^sin 

2 ( 27 ^~ 1 )( 27 ^ + 3 ) 


Um den Einfluß der Abplattung beurteilen zu können, soll nun die Schwere auf der Oberfläche 
einer Kugel gleicher mittlerer Dichte und gleichen Volumens eingeführt weiden, so zwar, daß 

G„ — — k^as 


und 

ist, dann erhält man 

gn ^^0 _ ^0 ^^0 


a = o\\ - —^^‘sin^s 
‘ 6 


r ^ 

äßd 277+1 5 


-sin- e 


'4 


1 


15 5 2(277 + 1) 


3 (2 77^—1) (277^+1) 


(2/7- 1) (277 + 3) 


1 


Im Falle in ~ 0, das heißt im Falle bloßer zonaler Eigenschwingungen ist dieser Einfluß leicht zu 
beurteilen. Die Gleichung 

Bn^Vgn — 0 


lautet dann 




4to2 


1 - 


2/7 (//+ 1) 
(277 - 1)(2^7 + 3) 


■ sin^ s 


1 


-77 (77+1) 


277-1-1 

GJi. 


77-—1 77(77 + 2 ) 


2/7-1 2/7 + 3 

3 Sr. 


1 - 


2/7+1 5 


4co^a^ 

Bezeichnet man mit X^ den Wert für s = 0, so ist 
1 


2/7 ( 77 +1) 


15 (277~ 1)(277+ 1)(277 + 3) 5 


—- sin^ £ 


aq 

4co"- 


und 


4 CO*' 


4(0' 


2 77 +1 
2/7 ( 77 + 1) 


7 / 2—1 77(77 + 2) 

2/7—1 2/7 + 3 


+ 77 ( 77 +1) 

4 03'- ä- 


(277- 1) (2/7 + 3) 4iü- 


^0 


1_ - .AL 

2/7+1 5 


ln ( 77 + 1) 


4(0^2 \15 (277— 1) (277+ 1) (2/7 + 3) '5 


sin- 











































Dynamische Theorie der Gezeiten, 


427 


Der rechts stehende Ausdruck ist aber jedenfalls positiv; negative Größen kommen nur in den von 
abhängigen Gliedern vor; in diesen hängt das Vorzeichen der von 5 unabhängigen Teile von 

2n(it+l) 4 

(27^-1)(2«+3) 15 

8 7 

ab, eine Größe die für n = 1 bis = oo Werte zwischen H-und H-annimmt, während der 

15 30 

s 4 s 

Koeffizient von zwischen 0 und — liegt. Da — gleich 0'18093 gesetzt werden kann, so ist die obige 
5 15 5 

Behauptung unmittelbar ersichtlich: X wird durch s vergrößert, woraus folgt: 

Die Abplattung verkleinert die Periode der freien zonalen Schwingungen. 

Was die numerischen Werte dieser Perioden anbelangt, so wurden diese von S. S. Hough für s = 0 
ermittelt lind zwar für verschiedene Tiefen einer die Erde gleichmäßig bedeckenden Wasserschichte. Es 
sollen hier nur die Ergebnisse angeführt werden, die sich auf jene beiden Annahmen beziehen, innerhalb 
deren die mittlere Tiefe der Wasserschichte liegen muß, die der Gesamtmasse der Ozeane bei völlig 
nivelliertem Erdkörper entspricht. 

Setzt man 

4w 40 ' 20 ' 

was für Jiq den Werten von 2200;;7, beziehungsweise 44007/7 entsprechen würde, so ergeben sich folgende 
Perioden, ausgedrückt in Stunden und Minuten (Sternzeit). 

nt—O Perioden 


n 

für/?o=2 

• 2 km 

— 

4 • 4 hn 

1 

30'' 

29’-3 

25'' 

p 

CO 

2 

18 

3-5 

15 

11-0 

3 

14 

15-2 

11 

54-3 

4 

12 

13-5 

10 

0-7 

5 

10 

50-0 

8 

38-3 

6 

9 

43-5 

7 

33-8 

7 

8 

47-5- 

6 

42-3 

8 

7 

59-6 

6 

0-5 

9 

7 

18-6 

5 

26-1 

10 

6 

43-3 

4 

57-4 

11 

6 

12-9 

4 

33-2 

12 

5 

46-4 

4 

12-5 


Zur Abschätzung des Einflusses 
mittlere Tiefe mitgeteilt werden: Für 

n = 

von s sollen noch folgende Daten für die erste angenommene 

. , AB , xg 

ist- und- 

4(0^ 4co^ 410^ 

1 

0-15491 

4-0-111s2 

2 

0-44155 

+ 0-215s2 

3 

0-70890 

+0-321 s2 

4 

0-96357 

+0-372s’“ 

5 

Wo X(, s “ 0 entspricht. 

1-2270 

+0 -436 52 
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Für den relativen Betrag - erhält man: 

0*364s2, 0-244s2, 0*216£S 0•193£^ 0*177£S 

er konvergiert, wie aus den obigen Formeln hervorgeht, gegen die Größe 

0 * 117 £ 2 . 

Die Fälle, für welche ist, die sich also auf die auch sektorial geteilten Eigenschwingungen 

beziehen, lassen sich hinsichtlich des Einflusses der Abplattung auf die Schwingungsperioden nicht so 
einfach behandeln, da der Mechanismus der Auflösung der Bedingungsgleichung ein ungleich kom¬ 
plizierterer ist. Es sollen auch hier nur die Korrektionsglicder für dinige der wichtigeren Spezialfälle 
angeführt werden. 

S. S. Hoiigh hat für verschiedene mittlere Tiefen die Perioden einer Gruppe dieser Schwingungen 
berechnet. Für die Tiefe 2*2 hUy die der am nächsten kommt, die dem Ozean als gleichmäßige blüssig- 
keitschichte entspricht, wurden von ihm die folgenden Zahlen gefunden: 


11 

in 

4 0)2 

oo 

Periode 


1 

1 

0-1024 

+ 3-125 

id 



2 

1 

0-6673 

+ 1-224 


14 

41“ 


1 

0-2418 

- 2-034 


24 

24 


2 

0-4454 . 

+ 1-498 


17 

59 1 


2 

0-0968 

- 3-215 


38 

34 

3 

1 

0-0010 

+ 31-995 

16 

0 



2 

0-0070 

+ 11-969 

6 

0 


4 

1 

1-0697 

+ 0-967 


11 

36 


1 

0-8312 

- 1-097 


13 

10 


2 

0-9866 

+ 1-007 


12 

5 


2 

0-6885 

- 1-205 


14 

28 1 

5 

1 

0-0004 

+ 52-854 

26 

10 



2 

0-0002 

+ 24-456 

12 

5 

! 

i 

Es wurde nun aus 

der Gleichung 




1 

1 




Bn+Sn + 0 







die Variation öo von und SP von der Periode, die aus s = 0 folgt, ermittelt für mäßige Beträge von L 
deren zweite Potenzen als zu vernachlässigende Größen vorausgesetzt wurden, und zwar für die ’ 
leie hen Argumente wie die obigen Resultate. fl 

Die erste Kolumne gibt die Koeffizienten von s- für So und die zweite die entsprechenden Faktoren | 
von 0 P, ausgedrückt in Sternzeitminuten: 


n 

m 

Sa; s2 

IP\Z^ 

1 

1 

-1-434 

-1033 

2 

1 

-0-176 

- 127 


1 

+0-954 

- 687 


2 

-0-293 

- 211 


2 

-0-231 

+ 166 
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i 


i 


n 

m 

8 0; 

8 P: £2 

3 

1 

+ 5-181 

+ 3731 


2 

+0-326 

+ 235 

4 

1 

-0-267 

- 192 


1 

+0-338 

- 243 


2 

+0-127 

+ 91 


2 

+0-234 

- 168 

5 

1 

+ 1-914 

+ 1378 


2 

+0-729 

+ 525 


Sind auf diese Art die Perioden ermittelt, so können für jeden Wert derselben die Größen kn und hn 
mit Hilfe der obigen Kettenbruchentwicklung angegeben werden und damit auch die Verhältnisse der 
unbekannten Koeffizienten An, so daß also für jede Schwingungsart, die durch X und m charakterisiert ist, 
die Amplituden für die einzelnen Breiten durch eine einzige willkürliche Konstante ausgedrückt erscheinen. 

Bezüglich der numerischen Ergebnisse sei auf die oben zitierten Arbeiten S. S. Hough’s verwiesen, 
die auch in Poincare’s »Theorie des marees« (Lecons de Mecanique celeste, t. III) zu finden sind. 


8. Erzwungene Schwingungen. 

Da hier X eine gegebene Größe ist, so sind die Koeffizienten a, b, c in dem Gleichungssystem der A 
ebenfalls von vornherein bestimmte Größen und es lassen sich daraus die An unter Berücksichtigung der 
Bedingung 

lim An ~ 0 

n = oo 


eindeutig bestimmen. 

Nun ergibt die Entwicklung der Störungsfunktion ein Aggregat von Gliedern der Form 

C„„n Pn,m (cOS &) e‘'"f. 

und zwar derart, daß jedem X nur ein einziges Glied entspricht, zu jedem X also nur ein Wertepaar nundm 
gehört. Es sei dieses 5 und in. Dann sind offenbar sämtliche An, für welche n nicht gleichartig mit 5 ist. 
Null, das Gleichungssystem der An für welche die n gleichartig mit 5 sind, lautet — je nachdem m gleich¬ 
artig mit .9 ist oder nicht: 

km Am G;i-(-2 Ani-\-2 Odcr. hjny\ Am-\-\ 4“ Cniy2> Ani-\-^ 0 

as—4: As—^ + bs~2 As—2-^^s — 0 

as-.2As—2'^bsAs-\“Csy2Asy2 — Cs 

Us Ais 4- ^s-(-2-^54-24-^5 _f_4 A.s_|_4 = 0 
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P'uhrt man wieder statt der A als Unbekannte die Verhältnisse /z„ und kn ein, so wird das System 
der Bestimmungsgleichungen für diese dasselbe sein wie für die freien Schwingungen mit Ausnahme der 
Gleichung, die hs enthält, demnach: 

hs—4i 4" ^5-2 +^5 — ö 

G 

hs~2 + ^5 + ^ 5 - i _2 —- 

Äs 

hsA- bs-\-2^ks-^4 = 0 


Die früher für die h und k gefundene Kettenbruchentwicklung wird also auch hier noch gelten, so¬ 
lange in dem Eliminationsprozeß die Gleichung mit der von Null verschiedenen rechten Seite nicht ent¬ 
halten ist, also zunächst für die Werte h bis hs-A und für die k von ^ 5+4 an, dann ergeben sich ks und hs 
aus der ersten und dritten der oben angeschriebenen Gleichung und aus der Indentität 


7 ^5—2 ^5 J 7 ^S +2 

hs -^2 — -und ^ 5_|_2 =- 

ks hs 

die linke Seite der mittleren Gleichung, so daß mit diesen Werten nun 

Q 


A=i- 


gefunden wird. 
Dann ist 


hs—2 4- ^5+^ 5+2 


As—^ — 


hs- 


s -2 


2 


^5—4 — 


hs —2 hs~A 

—2 ^^5—4 


^5+2 

<^5+2 G+4 


WO also sämtliche hier auftretende h und k aus der angegebenen Kettenbruchentwicklung erhalten werden. 
Man hat demnach folgendes Resultat: 

Enthält die Störungsfunktion ein Glied 

Cn,in Pn.m(cOS d) 
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so gibt dieses Anlaß zu Oszillationen, die definiert sind durch 



{ 


hji — 4, in hn 2, in 
—4, in ^11 — 2, in 


Pn~4,in (cos 4- 


(Ä»-2, m »< + ^« + 2, ,„)\/l — sin^ $ sin^ & 



Was das Größenverhältnis der einzelnen Glieder anbelangt, so kann zunächst bemerkt werden, daß 
das der Störungsfunktion entsprechende Glied mit Pn durchaus keine quantitativ ausgezeichnete Rolle zu 
spielen braucht: die Faktoren sind Funktionen von X, eine Größe, die ja bei diesen Schwingungen in 
keinem Zusammenhang mit dem Rang der Kugelfunktion steht. Es wird sich vielmehr in dem quantitativen 
Verhältnis ein anderes Phänomen bemerkbar machen können, nämlich das der Resonanz. 

Ist das X der Störungsfunktion eine Größe, die sehr nahe einem X der freien Schwingungen kommt, 
so heißt das nach dem vorausgehenden, daß es ein n gibt, für welches die Gleichung 


eine Wurzel in der Nähe von X besitzt, daß also für dieses X einen sehr kleinen Wert annimmt, daher 
auch das entsprechende An, wie aus dem System der Bedingungsgleichungen folgt, dementsprechend die 
übrigen Koeffizienten quantitativ übertreffen muß. 

Beschränkt man sich bei der Entwicklung der Störungsfunktion auf jene Glieder, die allein einen 
merklichen Beitrag zu dem Phänomen der Gezeiten liefern können, so kommt nur der Fall = 2 in 
Betracht und es handelt sich nun für die quantitative Beurteilung der einzelnen Oszillationen darum, ob 
die zugehörigen X der Störungsfunktion in der Nähe einer der Wurzeln 



9. Die Störungsfunktion. 


Ist in\ die störende Masse und A ihre Distanz von einem Flüssigkeitselement, so ist 


k'^ m ^ 
~ A 


Die den Koordinaten p, {)■, (p des letzteren entsprechenden elliptischen Koordinaten von itu seien 

9v fi 


und 


c 


~ sin s. 


dann ist nach bekannten Entwicklungen 


oo 



0 


11 



(/cots).QW(^cots.) 


Pj;;) (cos 0-) (cos cos in (cp^ — cp). 


(?;4-;/;)! (n — ni )! 


wo 
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Führt man wieder die Reihen und und die Neumann’schen Zugeordneten ein, so wird 


n 

T-oo = 2 yi . _ 

C ^ (ll+ 111)1 


Dabei ist also 


i?S«)(s).S(“)(s,)=r. 


cos” 


?” _ 

p” + ‘ COS’' + '”+' Sj 


1 + 


(n — tu) {ti — iii — 1) 
2.{2n~\) 


tg2 3+ . . . . 


(//-h///4“ 1) (?/4* ?//4“ 2) 

1 - ^-~tg2 £^+ . . . . 


2.(2//4-3) 


Da die Glieder u — O und ii — \ für die relative Beschleunigung keinen Beitrag liefern, — aber als 


kleine Größe betrachtet werden kann, von der nur noch die zweite Potenz berücksichtigt zu werden 
braucht, so kann 

s = r- r(2) 


gesetzt werden. Dann ist aber auch — zu vernachlässigen und wegen 

pi' 


= Pi i /1 --^cosä a-^ 


\ 


PI 


tgo = cot / \ — £L , 
\ PI 


wo r die geozentrische Distanz und 8 die Deklination von bedeutet, diese Größen für und 90°—0-^ 
zu substituieren. Dann wird 


S 




—- (3 — 2 sin^ s) fcos‘^ {)•-—] fsin^ o 

I ^]\ . 



4 - 2 cos s sin '8- cos 8* sin 8 cos 3 cos (rp — 
4 - -- sin^ 8* cos- 8 cos 2 C-f— • 


?~?i offenbar der Stundenwinkel des störenden Körpers für den Meridian Ist / die Sternzeit im 
Nullmeridian, von dem aus die p gezählt werden, so ist 


p — pj^ := co/ + p — a, 

wo a die Rektaszension des störenden Körpers ist, 

•Es handelt sich nun zunächst darum, die Variabein r, a und 8 als P\inktionen der Zeit und der Bahn¬ 
elemente des störenden Körpers auszudrücken, von dem man annehmen kann, daß er eine Keppler’sche 
Bewegung um den betrachteten Himmelskörper ausführt. 

Seien und Knotenlänge und Neigung seiner Bahnebene und co^ das Argument des Perigäums, 
bezogen auf das Äquatorsystem, und v die wahre Anomalie, so ist 

cos (a—cos 8 = cos (w^4-r) 
sin (a —cos 8 = sin (cOj4-r) cos 

sin 8 = sin (ü)^4-f) sin . 

Sind cfl, / und w die analogen Größen bezüglich der Ekliptik, der Bogen der Bahn zwischen den 
beiden Fundamentalkreisen und e die Schiefe der Ekliptik, so kann man nach den bekannten Trans- 
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formationsformein von einem System auf dcis andere übergehen, wobei von i nur erste Potenzen mit¬ 
genommen werden sollen, so daß 

sin /j cos = sin e-y i cos e cos £ 
sin sin = i sin ^ 

cos = cos e~ i sin e cos ß 
sin cosN= i cos ^+sin e cos ß 
sin i^ sin N — sin e sin 
= oj-hA^. 

Man wird sich bei dem vorliegenden Problem ziemlich weitgehende A'^ereinfachungen erlauben 
können, insbesondere sollen nur erste Potenzen der Bahnexzentrizität berücksichtigt werden. (Eine 
Verwechslung mit dem obigen das von dem elliptischen Bezugsystem herrührte und sofort gleich Null 
gesetzt werden könnte, ist ja nicht zu befürchten.) 

Ist die mittlere Entfernung und M die mittlere Anomalie, so kann darnach 

(r 

= 1-P3s, cos M 

\ / 

und 

V — il/+2 Sj sin M 

gesetzt werden. 

Nun sollen auch diese ersten Potenzen von i und nur in den von coi^ unabhängigen, also relativ 
langperiodischen Gliedern beibehalten werden, dann erhält man schließlich die folgenden drei Bestandteile 
der Störungsfunktion 

S ^ + So, 

wobei — unter tt die Länge des Perigäums und unter l die m'ttlere Länge in der Bahn verstanden — ; 

S. =- ^ -- (3 — 2 sin^ 3 ) { cos^ D---^ sin- e+ — sin^ e cos 2 / — 

4 ' rp ' i 3/132 2 

i 

-- [sin 2 e cos i^+sin 2 e cos (2 /—ß)] + 


+ 


3 1 7 

(1-sin- e) cos — - -sin^ e cos (/ + t:) H-sin 2 e cos (3/ - tu) 

2 4 . 4 


S. = 


3 

— k^ni^ — cos £ sin {> cos 9* sin e [sin (ooZ+'f) —sin (co^+rp — 2/)] 


/o 


(dabei ist sin^^ als von der Ordnung i angenommen). 




3 (ß“ 

— sin^ 9’ [sin^ e cos 2 (()o/ + cp) + 2 cos 2 (co/~hrp —/)] 

8 rl 


mit der gleichen Vernachlässigung. 

Wie man sieht, tritt in jedem dieser Bestandteile eine einzige Kugelfunktion auf, und zwar in Sm 
nurP2,m. 

Für 77/:::i;0 erhält man die Gezeiten langer Periode, das Hauptglied hat eine Periode die gleich ist 
der halben tropischen Umlaufszeit des störenden Körpers, involviert also etwa vierzehntägige Mond- oder 
halbjährige SonnenOuten. 

und geben Anlaß zu kurzperiodischen Gezeiten: enthält Perioden, die gleich einem Stern¬ 

tag oder um einen geringen Betrag länger sind, S., halbtägige oder nahezu halbtägige Perioden. 

DsnUschriflen der mathem.-iiaturw. Kl. LXXXlK. Bd. »^7 
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IO. Langperiodische Gezeiten. 


Hier ist m~0 und —eine kleine Größe. 

2 tü a 


Die Betrachtung der Koeffizienten 


Un 


bn 


Cji IIq 


4 ü)2 p- 4 to^ p- 4 0)'^ p- 

ergibt dann, daß diese so gut wie unabhängig von X sind, wegen der Kleinheit der von X abhängigen 
Glieder. Die numerische Rechnung zeigt tatsächlich, daß die Resultate für X =*0 nicht wesentlich von 
jenen äbweichen, die man für den größten hier in Frage kommenden Wert von X, den der vierzehntägigen 
Mondflut, erhält. Aus demselben Grunde sind diese Größen auch merklich unabhängig von s, denn auch 
in gn kommt die Exzentrizität nur mit einem kleinen Faktor versehen vor, der gegenüber dem expliziten 
Vorkommen von s in der Störungsfunktion als belanglos erscheint. 

Setzt man wieder 


so findet man für 




hQ = 2-2 km 

Go K 



4(1)2 

log di- 

A 

=: 9„7979 

, A 
log-— = 

A 

9„4Ö01 

logd^ 

= 9„033 

log-^ = 

8 „88. 


40 ’ 


log -d8_ — 9„22ll 


Die Höhe C der Gezeitenwelle wird also für diese langperiodischen Bewegungen gegeben sein durch 
eine Reihe nach Pn,o (cos ^), deren erstes Glied uz=z2 das quantitativ bedeutendste ist. 

Es ist also 


C \/1 — sin^ £ sin^ff = ^2 + ^4 ^^4+.) 


wo nach dem früheren 




K+c, —^ • 


Man findet gemäß den angegebenen Voraussetzungen 


A, — 19-28. 


C, 

4ü) 2p-' ■ 


wo 


c = -- Hm, ^ 


1— sin-£|sm“^, 

3 ' 


A., - 5-423. 


/;q 

(i)” / 


l 3 


sin- eil -sin^s 


so daß 
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Was zunächst den Einfluß der Abplattung anbelangt, so sieht man, daß die einzelnen Amplituden 
proportional dem Faktor 

1 — sin‘^ £ f — “ sin^ 

\ 3 2 

verkleinert werden, also im Maximum an den Polen und am wenigsten am Äquator. 

Für die halbjährige Sonnenflut hat man 

in^ uz 1 und Tq = 1 

zu setzen und man erhält, wenn für 

d’ = 23^^ 2f 8" 


und der Annahme gemäß h^^ i= 2200 m gesetzt wird, 


yjg = 0*0142 m I 1 —^ sin^s^ . 


Für die vierzehntägige Mondflut hat man 


k'^ m. 


Ul 


wo [X die Mondmasse in Einheiten der Erdmasse und die tropische Umlaufszeit des Mondes in Stern¬ 
tagen bedeutet. Für 

1 


findet man 


demnach 


[ 11 =- und 27*3291 

81*45 


. i ‘2 

— 0*0391 m j 1-sin‘-^£ 


= 2-76 

A, (O) 


Für größere Werte von nehmen die Amplituden rascher zu, weil der Faktor von h^ in selbst 
mit Hq zunimmt; so wird er etwa für die vierfa^jhe der hier angenommenen Tiefe 2* 135 mal größer und 
die Amplitude der Mondfluten für Iiq =8*8 km wird demnach sein 

A, = 0*2654 m |^1 
\ 



II. Ganztägige Gezeiten. 

Die von herrührenden Glieder in C haben Perioden, die gleich oder nahezu gleich einem Stern¬ 
tag sind. 

Es kann gezeigt werden, daß die Flutbewegung dieser Art unter allen Umständen nur außerordentlich 
kleine Amplituden haben kann. 

Der hier vorliegende Fall ist in allgemeinerer Fassung der folgende. 

Die Entwicklung der Störungsfunktion habe für ein bestimmtes m und X nur ein einziges Glied, und 
zwar soll für dieses n =7;z + 1 sein, so daß also 

S = Pin (cos ö*) e^”‘ ^ e^^"^ 
ist. 
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Dann werden sämtliche sein und die A,n-\- 2 k^i werden dem System in Nummer 6. zu 

genügen haben 

kjH \ 4 “ hn-{-2^ 

-4;;/+3‘4-^/;2+4 — 6 

A^7;2 m-\-l Ajn-l^] ~ 0 

k}n-\-\ iY;; 2^2 Djii-^2~^hn-\-3Aiii-\^Z — 0 


Nun sei X so beschaffen, daß A'^ =0 wird, das heißt 


oder 


in (nt-}-—am—s\n^ = 0 

. 2 Cf) 

nt cos‘^ s + 1 ’ 


dann ist N[,„^i z=0 und dem obigen Gleichungss 3 ^stem wird offenbar genügt durch 


Du 


IZjj 


j Dju^2 — — Djji^4. — u4;,2_j_5 — .... — 0, 


das heißt, das betreffende Glied der Störungsfunktion gibt zu keiner Gezeitenbewegung Anlaß. 

Im vorliegenden Falle ist nt— 1, n = 2 und X = co oder sehr wenig davon \^erschieden. 

Vernachlässigt man s. so wird die Bedingung —erste Glied von streng, für das 
zweite mit großer Annäherung erfüllt sein, und man erhält, da dann h= ist, das bekannte Theorem von 
Laplace, demzufolge bei einem Ozean von konstanter Tiefe keine ganztägigen Gezeiten existieren. 

Bei kugelförmiger Erde kann das für die Sonnenflut ohneweiters angenommen werden; mit geringerer 
Annäherung gilt es für die Mondflut, da hier der Faktor von / in 2/ nicht mehr unmerklich klein gegen 
CO ist. 

Nimmt man aber Rücksicht auf die Erdabplattung,^so ist das auch für X=:co nicht streng richtig. 

Die quantitativen Verhältnisse ergeben sich dann in folgender Weise. 

Angenommen es sei X=: co + wo t eine kleine Größe bedeutet, so ist, wenn zunächst s = 0 voraus¬ 
gesetzt wird, 

. 

0 ) 

Das Gleichungssystem der A ist hier 

D ^ ■— il/,^2 4-/3 D .^ ~ 
k,D,-M,A,+I,D, = 0 


-A\D,+J,A, = 0 

/v% A2 — A3 A^ — 0 
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woraus folgt 


. ^ z D. 

— 2 -, 

to L 


worin Z)i = — gesetzt werden kann, da der wahre Wert von von diesem nur um Größen von der 

K 

z 

Ordnung — abweicht. 

0) 


Es ist demnach 


und, da hier 


ist, so findet sich 


4 = 




jy 

Kh 


3 T 2 

— 

2 




rl 


z k- in, h 

lo- — 

(0 rl co^ 


sin e 


sin e 


Das Glied in das ein von Null verschiedenes t hat, ist das mit dem Argumente oj/ + rp — 2/; 
T ist also negativ, die Amplitude dieser Gezeiten wird der des entsprechenden Störungsgliedes entgegen¬ 
gesetzt sein. (Die Koeffizienten . . . nehmen sehr rascli ab.) 

Da für die Sonne — =-^-> für den Mond hingegen = — ^ ^~ , so folgt für das Verhältnis 

(0 182-62 13-66 


(A.^^:(A,)q^J^\ = 29-33, 

\^C / 

so daß die ganztägigen Sonnenfiuten tatsächlich unmerklich sind. 

Die Höhe der Gezeiten nimmt auch hier mit der mittleren Meerestiefe zu, aber bei exakter Rechnung 
auch wieder rascher als nach bloßer Proportionalität. Für sehr große Tiefen tritt übrigens hier der oben 
erwähnte Fall der Resonanz ein, das heißt, die hier bestehende Periode der erzwungenen Schwingung 
kommt mit zunehmender Tiefe der Wurzel einer Periodengleichung der freien Schwingungen immer 
näher. Die Periode der hier auftretenden Mondflut beträgt ja ungefähr 1^ Tü für u — l und 1 erhielt 
man bei der Tiefe 2'2 hu für die freie Schwingung eine Periode von 1^ 13^G Annahme größerer Tiefe 
wird die Periode kürzer, erreicht bei achtfachem Betrag der mittleren Tiefe den Wert 1^ 3’^ und nähert 
sich immer mehr der Mondflutperiode. Die Rechnung ergibt, daß für diese Tiefe die Amplitude der ganz¬ 
tägigen Mondflut das Neunzehnfaehe der für h^ — 2'2hn beträgt. 

Was nun den Einfluß der Erdabplattung anbelangt, s6 folgt zunächst 



CO 


und für die Amplitude. 

l2z \l 1 \ 

A 2 = -sin“^ sj l^l-h sin- sj (AA.=o 

woraus einerseits hervorgeht, daß mit Rücksicht auf das Vorzeichen von z diese Gezeiten vergrößert 
werden, anderseits daß auch die Gezeiten mit genauer Sterntagperiode einen von Null verschiedenen 
allerdings sehr kleinen Betrag erhalten, der etwa von der Größenordnung der Sonnenfluten beiü<ugel- 
förmiger Erde ist. 

Es soll übrigens hier noch bemerkt werden, daß die beiden bis jetzt untersuchten Gezeiten, die 
langperiodischen und die ganztägigen wesentlich durch die Abweichung der Bahnebene des störenden 
Körpers vom Rotationsäquator bedingt sind und nicht existieren würden, wenn die betreffende Bahn in 
der Äquatorebene liegen w'ürde. 
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12. Halbtägige Gezeiten. 


Die von S.^ herrührenden halbtägigen oder nahezu halbtägigen Gezeiten bilden die quantitativ weit¬ 
aus überwiegenden Komponenten des ganzen Phänomens und verleihen diesem den auffälligen charak¬ 


teristischen Verlauf. Es sind dies jene Erscheinungen der Ebbe und Flut, die schon aus der statischen 


Betrachtungsweise unmittelbar erfolgen — wenigstens ihrem typischen Verhalten nach — und die auch 
der oberflächlichen Beobachtung ohneweiters zugänglich sind. 

Um eine Vorstellung von der Größenordnung dieser Gezeiten zu erhalten, soll zunächst wieder die 
Abplattung vernachlässigt werden. 

Es soll hier der gebräuchliche Vorgang eingehalten werden, die Amplituden im Verhältnis zu denen 
der statischen Gezeiten anzugeben, die erhalten werden, wenn in den Fundamentalgleichung die als Faktoren 
auftretenden X 0 gesetzt werden. Sind diese Gezeiten gegeben durch 


C\/1 — sin^ s siiV-^ d- = 1 3(„, ^ (cos 0-) , 


so reduziert sich die aus der Oberflächenbedingung resultierende Bestimmungsgleichung der Koeffizienten 
auf 


gn — Cn 


es entspricht also jedem Glied der Störungsfunktion nur ein einziges Glied in C- 

Das in auftretende Hauptglied, mit dem Argumente 2 (w / + ^ — /), gibt, weil 


P 2,2 (cos = 3 sin‘^ ^ 


für 



demnach 



16 



oder für s = 0 



wo — = 0* 18093 ist. 


•^1 


Nimmt man wieder Jiq = 2'2km demnach 


Uo 1 

4oV ~ 40 ’ 


so wird 


0*022286 %, 

Die dynamische Theorie der Gezeiten hat ergeben 
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Die oben in 7. angegebene Kettenbruchentwicklung für die Größen 


-(-2 — + 2 


^7^ + 2 
^77 


gibt für n—2, m = 2 


log 


4(02f/^ 


~ 8-71588 


ist unmittelbar der entsprechenden Formel zu entnehmen, man erhält schließlich 


_l-9476. 

9f,, 

Um das tatsächliche Amplitudenverhältnis zu finden, hat man in derselben Weise die Reihe der 
Koeffizienten Aq, . , , . z\x bilden. Man findet 

log A =0-03810 

A 

log ^ = 9„5826Ö 

log ^ 9„22417 

^6 


und daraus 

A, - 1 • 9476 
A^=: -2-12617 %, 
Aß =1 +0-81331 % 
A^ = -0-13628 21.. 


Die Amplitude der statischen Gezeiten ist 


die der dynamischen 


wo 


“^2 ^2 >'2 

[A. P .,2 (cos {)’) + Aj^ (cos \)')+ . . , • ] ? 


P„,2(COS^) 


1.3....(2n-l) 

{ 11 - 2)1 


sin^ ^ 


cos^^~2 ^ __ 


(n—2) (n — 3) „ 40 . 

2.(2^^-l) 


+ 


(n-2) . . .(n-5) 

- -^^cos^'-<5 

2.4.(27/-l)(2n-3) 


Die numerische Ausführung ergibt, daß für den Äquator, 'O-mOO'', die Amplitude der dynamischen 
Gezeiten sehr nahe das Dreifache von ist. 

Um absolute Beträge zu erhalten, hat man in c die entsprechenden Elemente der störenden Körper 
einzuführen. 
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Es ist 

log c '0 " 3‘(37070 - 10 
- log C(r — 4-01173-10 

woraus für den angenommenen Wert //,) folgt 

(O) = 0-04625 m 
% (C) —0-10142 m, 

so daß also am Äquator der Maximalbetrag der halbtägigen Gezeiten 

0-886 m 

betragen würde, bei einem Ozean, der den ganzen Erdkörper mit einer gleichmäßigen Tiefe von 2'2 km 
bedeckt. 

Bei der Durchführung dieser Rechnung wurde der durch l auftretende P'aktor von t gegen co als 
belanglos vernachlässigt. 

Das gleiche gilt auch von dem Einfluß der Abplattung für die Verhältnisse auf unserer Erde. 

Will man ganz allgemein diesen Einfluß bei elliptischen Rotationskörpern berücksichtigen, so sieht 
man zunächst, daß der betreffende Teil der ätörungsfunktion unabhängig von s ist, die A also nur ver¬ 
möge der Änderung der Größen Z?«, Cn zu variieren sind. 

Über den Sinn dieser Variation läßt sich von vornherein nichts aussagen, es hängt das ganz von 
den speziellen Voraussetzungen des Problems ab. 

Setzt man ähnliche Verhältnisse wie auf der Erde voraus, so kann die folgende Betrachtung darüber 
Aufschluß geben. 

Es soll sich in erster Linie um die Variation der beiden quantitativ bedeutendsten Koffizienten 
und handeln und vorausgesetzt werden, daß und der Faktor von in Aq so klein ist, daß Aq nicht 
mehr merklich von s beeinflußt wird. Dann hat man zur Bestimmung von A., und A^ 

h.) ^2 A^ = C., 

Li.^ A.2 H“ A^ — Cq A^ 

und, da die beiden rechten Seiten voi\s unabhängig sind, so hat man für die gesuchten Variationen 6 A 

a., ö A.y4- c^o A^ — — s- A., f ß., -i- Y 4 j 

‘ \ ’ A.J 

a.yZA.y4-b^oA4^— — A^ -i-ß^ ^j 

wo die oc, ß, Y die Koeffizienten von sin^s in den entsprechenden Größen a, b, c sind. 

a und c sind wesentlich positive, a und y negative Größen, die mit wachsendem ;; und konstantem ni 
(— 2 ) gegen Null konvergieren; in | 

^n—Sn4-Bii \ 

ist wesentlich positiv und konvergiert gegen eine endliche Grenze, negativ und gegen Null kon- « 
vergierend, es ist also hji positiv oder geht einmal von einer negativen Anfangsreihe zu positiven Werten " 

über, für ß besteht ein ähnliches Verhalten. Tritt dieser Wechsel schon zwischen n — 2 und ;/= 4 ein, so ” 

ist die Determinante wesentlich negativ und es wird sich das Vorzeichen vonoA^undoA^ leicht ent- j- 
scheiden lassen. 

Im Falle der Erde tritt das tatsächlich ein bei der oben angenommenen Tiefe der Wasserschicht. f 
Alan erhält für die Variationen oA Größen, die mit A^, und daher auch mit A^ verschiedenes \'or- : 
Zeichen besitzen, das heißt, die Amplituden werden kleiner. Es ergibt sich 

A, z= - 1 • 9476 (i - 4 * 43 s'^) 

A^=:— 2-1262(1 - 5 96 3-051,. 









